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O ■ ABSTRACT 

f^*) \ GMOS optical long-slit spectroscopy at the Gemini-North telescope was used to clas- 

t— I • sify targets from the Supernova Legacy Survey (SNLS) from July 2005 and May 2006 

- May 2008. During this time, 95 objects were observed. Where possible the objects' 
• i-H , redshifts (z) were measured from narrow emission or absorption features in the host 

■ galaxy spectrum, otherwise they were measured from the broader supernova features. 

We present spectra of 68 confirmed or probable SNe la from SNLS with redshifts in the 
range 0.17 ^ z ^ 1.02. In combination with earlier SNLS Gemini and VLT spectra, 
we used these new observations to measure pseudo-equivalent widths (EWs) of three 
spectral features - Can H&K, Sin and Mgn - in 144 objects and compared them to 
the EWs of low-redshift SNe la from a sample drawn from the literature. No signs 
of evolution are seen for the Can H&K and Mgn features. Systematically lower EW 
Sin is seen at high redshift, but this can be explained by a change in demographics 
of the SNe la population within a two-component model combined with an observed 
correlation between EW Si II and photometric lightcurve stretch. 

Key words: supernovae - cosmology: observations 
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1 INTRODUCTION 



Since the discovery of cosmic acceleration (e.g. iRiess et al.l 
1 19981 ; IPerlmutter et all [l999), Type la supernovae have 
played a vital role as distance measures in the universe, up 
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to and beyond z = 1 th rough surveys such as the Supernova 
Legacy Survey (SNLS Ustier et all l2006l l and ESSENCE 
jWood-Vasev et al.|[2007F 7 This is made possible by the fact 
that SNe la appea r to be " s tanda rdisable" candles, a prop- 
erty first noted by I Phillips] |l993) and la ter parameterised 



in a number of differen t ways (for example IRiess et al.l 
IPerlmutter et all 1 19971 ; IWang et all 120031 ; iGuv et al 



1996; 



2005; 
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IJha et al.ll2007h . However, the use of the same lightcurve 
calibration, e.g. stretch, a measure of the relative width of 
the lightcurve, and colour, at high and low redshift relies on 
the fact that there is no evolution in the properties of the 
supernovae with time that cannot be accounted for by the 
calibration itself. 

The properties of high- and low-z objects have been 
compared in a number of ways. One way is to compare the 
line profiles and eject ion velocities of features in the spec- 
tra. iHook et al.l (|2005l ) measured the ejection velocity from 
the blueshifted Call H&K feature in 11 high-z supernovae 
with z > 0.38 and compared this to low-z objects. They 
found that at high-z the ejection velocities and phase evolu- 
tion of the spectra are consistent with the low-z sample and 
concluded that there exists a sample at high-z w ith similar 
properties to SNe at low-z. iBlondin et alj (200a) measured 
the wavelength of maximum absorption and peak emission 
for the Can H&K, Sin 6150 and Sn "W" features and found 
them to be similar in high-z and low-z objects. They also 
noted evidence of a double absorption feature in the Can 
H&K region which is seen to occur occasionall y in low-z 
super novae and has been attributed to Sin (| Howell et al.l 
l2006h . 

Since the advent of large surveys producing high num- 
bers of spectra, it has been possible to use these to make 
composite spectra, which are representative of the average 
supernovae spectrum in certain redshift and phase bins. 
This has been done using Keck spectra of SNLS super- 
novae by lEllis et all (|2008l . h ereafter E08 ) with z = 0.5, 
and for ESSENCE spectra by iFolev et al.l (|2008l , hereafter 
F08) . The studies used different methods to account for host 
galaxy contamination, but neither found definitive evidence 
for an evolving supernova population. |E08| created a compos- 
ite in order to study possible effects of metallicity variations 
in progenitors in the rest-frame UV spectra. This is where 
the absorption is dominated by metal line-blanketing effects 
and the metallicity differences are expected to be the most 
prominent, although the theoretical predictions are not in 
agre emen t (|Hofiich et al.lll998l ; iLentl et al. I l2000l . for exam- 
ple). IE08| saw a large amount of variation in the spectrum 
in the UV region, more than is predicted and can be ac- 
counted for by dust, and also changes with phase that are 
not predicted. However, due to a lack of local UV spectra of 
supernovae they are unable to draw any firm conclusions as 
to evolution. 

|F08| do not perform galaxy subtraction, but instead con- 
taminate their low-z spectra with a known galaxy spectrum 
in an attempt to match the host galaxy contamination of 
their high-z sample. They find evidence of changes in the 
Fell spectral feature which may indicate lower 54 Fe produc- 
tion in a lower temperature environment brought on by a 
lower metallicity. They attribute this not to evolution, but 
to changing demographics as the universe at higher redshift 
was less metal-rich compared to the universe today; how- 
ever, again their conclusions are limited by the low-redshift 
sample of spectra. 

The composite spectra from lEOSl was included in 
ISullivan et al.l (|2009T ) which compared their medium redshift 
composite to a new low-z composite (z < 0.05) and a high- z 
composite (z > 0.9) based on data from lRiess et al.l (|2007h . 
They observe a clear decrease in the size of spectral features 
of intermediate mass elements as redshift increases, but this 



is a change which is consistent with a change in demograph- 
ics, rather than an evolution of SNe la. 

Other studies have focussed on measuring the pseudo- 
equivalent widths of spectral features in the high-redshift 
spectra and comparing these to trends with phase in the low- 
redshi ft sample. This method was first described bv lFolatellil 
(2004, hereafter F04) and used successfully on a small sam- 
ple of spectra to compare low- and high-redshift objects for 
5 distinct features in the spectrum. The blueshift of the 
absorption minimum in Call H&K was also measured and 
compared to nearby supernovae and no evidence for any sig- 
nific ant differences was fo und between the two samples (see 
also iGaravini et~al"1l2007h . 

An EW analysis of the first two a nd a half years o f Gem- 
ini SNLS spectra was published in iBronder et al. I (|2008l . 
hereafter B08) . Rest-frame EW measurements were made for 
the Can H&K, Sin 4000 (hereafter referred to as just Sin) 
and Mg II features as well as the Ca II H&K ejection velocity. 
The spectra of 55 high-z supernovae were compared to 167 
spectra from 24 low-z objects encompassing all three SNe la 
sub- types (normal, overluminous and underluminous). No 
significant difference was seen between the two samples for 
the equivalent widths of Ca II H&K and Si II or the ejection 
velocity. For the Mg II feature, a slight difference was noted, 
but it was not significant enough to confirm evolution and 
could just have been a stat istical effect. This is discussed 
further in Section 14.41 |B08| also found a relation between 
the EW Sill and Mb, the absolute B-band magnitude at 
maximum, which could be used to calibrate the supernovae 
in the same way as lightcurve stretch and colour. 

With the well-measured light curves of the SNLS su- 
pernovae it h as also been possibl e to compare photometric 
properties. In lConlev et alj l|2006l ). the rise-times of low- and 
high-z supernovae are compared and they find no compelling 
evidence that there is a difference between the two popula- 
tions. 

However, as was concluded by iFOSl . there is to be ex- 
pected some change in photometric properties due to chang- 
ing demographics. Due to the two-channel progenitor mode l 
ilScannapieco fe Bildstenl 120051 ; iMannucci et al.ll2005l . 120061 ; 
ISullivan et al.| [2006). the "prompt" component which are 
brighter supernovae with broader light curves should dom- 
inate at higher redshift, whereas at z = 0, the "delayed " 
component wi l l be more important. iHowell et alj |2007); 
ISullivan et al] (|2009h calculate that this should result in a 
6% increase in light curve width between z = and z = 1.5. 
Using observational data from ISullivan et alj (|2006l ) , they 
measure a percentage change in stretch of 8.1±2.7% between 
0.03 < z < 1.12 which is consistent with c hanging demo- 
graph ics. The spectral evolution observed in ISullivan et alj 
(2009) is consistent with this. 

In summary, no firm evidence of evolution in the spec- 
tral properties of supernovae has been observed, but the 
studies are being hampered by the lack of local spectra to 
compare with. In this paper we present an analysis of EWs of 
the Call H&K, Sill and Mgn spectra l feat ures. This study 
significantly expands on the work of iBOSl using new data 
from SNLS with a different method of accounting for the 
effects of host galaxy light, as well as an expanded sample 
of low-z SNe from the literature for comparison. 

The remainder of this paper is organised as follows. In 
section 2 we describe the new observations. In section 3 we 
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detail the methods and measurements applied to the high- 
z sample of SNe la. Section 4 describes the low-z sample 
drawn from the literature and used for comparison. Section 
5 details the results of the comparison study and an investi- 
gation into the possible cosmological uses for SNe la spectra. 
These are discussed further in section 6 and our conclusions 
are summarised in section 7. 



2 OBSERVATIONS 

The targets were selected from real-time SNLS imaging 
which took place at the Canada-France-Hawaii Telescope. 
Fainter objects were usually observed spectroscopically at 
Gemini, with brighter candidates at either VLT or Keck 
(see Perrett et al. (in prep.) for more details on the real- 
time pipeline). 



2.2 Data Reduction 

The data were reduced using standard IRAF0 and Gemini- 
specificlfl tasks to perform bias-subtraction, fiat-fielding, sky- 
subtraction and chip mosaicing. A wavelength solution was 
applied to the 2D spectrum. When combining the 2D frames, 
bad-pixel masks were used to take account of the "missing 
data" in the gaps between the CCDs. A ID spectrum was 
extracted from the 2D spectrum and corrected for telluric 
absorption features. It was generally not possible to extract 
the target light separately from the host due to the small an- 
gular separation of the objects on the sky. A sky spectrum 
was also extracted from the 2D frames, and an error spec- 
trum was estimated from this by assuming that the noise 
is dominated by sky emission. Where appropriate in some 
of the earlier data which were observed at only one central 
wavelength, a linear interpolation was used to remove the 
chip gaps and the error spectra set to a high value in these 
regions to deweight them. 



2.1 Observational Set-up 

The new spectra presented here were taken in queue mode 
on the Gemini-North telescope with ToO status^]. The 
observational set-ups and data reduction me thods used 
for S NLS-Gemini spectroscopy using GMQS jHook et alJ 



2004) have been described in detail previously (I Howell et al 



2005 ; [Bronder et al.ll2008l ). In summary, a 'classical' setup 



with central wavelength of 680nm was used for obser- 
vations of brighter targets which were expected to be 
at a lower redshift, whereas nod-and -shuffle mode (N&S, 
iGlazebrook fe Bland-Hawthorn! 1200 ll ) with a redder cen- 
tral wavelength of 720nm was used for fainter targets. The 
only change t o the methods of observ ationa l set-up since 
iHowell et all l|2005l . hereafter H05) and lB08l are that from 
March 2007, first N&S and then classical observations were 
taken at two wavelength settings in order to cover the gaps 
in the spectra caused by the spacing of the three GMOS 
CCDs. In these cases, for the R400 grating and with 2x2 
binning, the change in central wavelength used was 50 nm. 
Due to the unavailability of the R400 grating, three targets 
were observed with the R150 grating which has a lower dis- 
persion. For these objects, a change in central wavelength of 
100 nm was used. 

The position angle of the slit was primarily determined 
by aligning the target and the centre of the host galaxy in 
order to try to observe narrow emission/absorption features 
in the host spectrum. The second constraint was the avail- 
ability of guide stars in the field. Occasionally, less optimum 
position angles were used or the object was moved from the 
centre of the slit. All of these considerations mean that the 
observations were not taken at parallactic angle, but this is 
taken into account in the analysis that follows. 

The new targets observed as part of the Gemini SNLS 
spectroscopic follow-up programme are listed in Table [ATI 



1 Observing programmes: GN-2003B-Q-9, GS-2003B-Q-8, GN- 
2004A-Q-19, GS-2004A-Q-11, GN-2004B-Q-16, GS-2004B-Q-31, 
GN-2005A-Q-11, GS-2005A-Q-11, GN-2005B-Q-7, GS-2005B-Q- 
6, GN-2006A-Q-7, GN-2006B-Q-10, GN-2007A-Q-8, GN-2007B- 
Q-17, GN-2008A-Q-24. PI Hook 



2.3 Redshift Determination and Object 
Classification 

The redshifts of the supernova and its host are ideally mea- 
sured using narrow host features such as [Om], [On], H/3 
and Ha emission, and H&K absorption. This typically gives 
an error of z ± 0.001, but if many lines are present then the 
error can be reduced further. If no host features are present 
then the redshift measurement relies on the wider supernova 
features. These are less precise due to the higher velocities 
involved and typically measure z ± 0.01 depending on the 
quality of the spectrum. 

A high-z SN la is classified by identifying, typically, the 
Sin feature at A ~ 4000 A in the spectra. Other features 
used are Sin A6150 and Sn. To aid classification, all spec- 
tra are also run through the SUPERFIT program 
which performs a ^-minimisation to best-fit the observed 
spectra to a low-z supernova spectrum combined with host 
galaxy light. Objects are classified as SNIa (certain or highly 
probable la), SNIa? (probable la), SNII (SN Type II), SN 
(supernova of unknown type) or SN? (probable SN of un- 
known type). To go alongside this, spectra are assigned a 
confidence index (CI) from 0-5 with 5 being a definite la 
and definitely not. Objects of CI = are typically core- 
colla pse supernovae. The definitions of the CIs are given in 

IhoI 

The redshifts, their errors and CIs of all the objects 
observed in this study are given in Table [ 



2.4 Sub-Type Identification 

A number of SNe la were observed to have spectral fea- 
tures similar to kn own p eculiar objects from their best-fit 
SUPERFIT matches |H05T ). Identification of these objects is 
important so they can be flagged and removed from later 
studies. The objects which were identified as having unusual 
spectral features during this period of observation are given 
in Table □ 



2 http://iraf.noao.edu/ 

3 http://www.gemini.edu/ 
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Table 1. Unusual SNe la identified in this sample from 
their spectral features only. The stretch measurements are 
made using the SiFTO dConlev et alj l2008ri light curve fitter. 
The references for the comparati ve objects are SNLS-03 D3bb 
llHowell et al.l l200rj). SN2001ay iHowell fc Nugent! l2004h and 
SN1999aa jGaravini et al]|2004f ). 



Object 


Stretch 


Comments 


07D3af 


0.943 ± 0.017 


03D3bb-like 

C II feature at 4200A 


07D3ap 


0.955 ± 0.018 


SN 2001ay-like 
High- velocity Sill 


07D3cr 


0.872 ± 0.043 


SN 2001ay-like 
Broad Call 


07D3dj 


1.118 ± 0.018 


03D3bb- or SN 1999aa-like 
Little Ca II and probable C II 


08D2aa 


1.128 ± 0.016 


SN1999aa- or SN 2001ay-like 
C II feature at 4200A 


08D2ad 


1.155 ± 0.021 


03D3bb-like 

C II feature at 4200A 



There are several key features which are distinct to 
each of the compa rison spectra. SNLS-03D3bb (SN 2003fg 
I Howell et al.|[200g) was an overluminous SNe la whose pro- 
genitor system is debated, but it is thought to be of super- 
Chandrasekhar mass. The spectrum showed low ejection 
velocities, as well as an unburnt Cn feature at 4200A. 
SN 2001ay l|Howell fc Nugentll2004r i has a high stretch value 
(s = 1.6), but appeared spectroscopically normal around 
maximum light, e xcept for high-velocit y Call and Sill fea- 
tures. SN 1999aa |g aravim et al.ll2007f ) has been seen as a 
transitional object between normal SNe la and the overlum i- 
nous objects typified by SN 1991T l|Filippenko et al.lll992l ). 
showing spectral features of both types of object. SN 1999aa 
spectra at maximum contained smaller amounts of inter- 
mediate mass elements and larger amounts of iron-group 
elements, and appeared peculiar until maximum light after 
which it appeared to be spectroscopically normal. 



3 METHODS AND MEASUREMENTS 

For the use of SNe la as cosmological distance indicators, it is 
important to verify that the objects we observe in the high-z 
universe have similar properties to those at low z, or, if there 
is some change in the mean properties l)Howell et al ] |2007l . 
for example), that these differences can be calibrated out. 
All p r evious work on high- z SNe la spectra 

iBronder et al 



20051 ; iBlondin et all 120061 : lEllis et all I2OO 
20081 ; iFolev et al.l 120081 ) have found no definitive evidence 



for evolution beyond what is expected from known changes 
in demographics. 

In this paper out of the 95 objects observed, 68 were 
confirmed as Type la or probably Type la events. Pr evious 
Gemini SNL S samples have published a further 41 l|H05h 
and 46 llBOsh. 



3.1 Sample Selection 

For the study here, all SNLS super novae obs erved with Gem- 
ini, including objects published in lH05l and lB08l . with confi- 
dence indices of 4 or 5 were considered. A cut was then made 
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Figure 1. Redshift histogram for the Gemini and VLT samples. 
The dotted line represents the mean redshifts of the two samples. 



based on the availability of final lightcurves and measure- 
ments of conditions on the night of observation, specifically 
hour angle (HA) and mean seeing. The HA and seeing are 
required for light-loss correction described in Section [3.3.31 
In the case of observations spread over more than one night, 
the larger value of the seeing was used. 

A comparison sample of spectra drawn from the first 
three years of SNLS- VLT spectra analysed in real-time 
(Basa et al. in prep) was also used. T hese objects were 
obser ved using the FORS1 spectrograph (jAppenzeller et all 
1998) on the VLT. Again, these objects were required to 
have CI=4 or 5 and have final lightcurves and a measure- 
ment of seeing for the EW study. In this case, the measure- 
ments of the hour angle on the night of observation were 
not required as FORS1 employs an atmospheric dispersion 
correction (ADC) which minimises the amount of light lost 
when not observing at the parallactic angle. 

Removing objects classified as peculiar (see Table [TJ, 
the sample for the EW study contains 149 objects, 82 from 
Gemini and 67 from the VLT. A histogram of redshift is 
shown for the final Gemini and VLT samples in Figure [T] 
It shows that the Gemini sample peaks at a higher mean 
redshift, as expected since Gemini was used to observe the 
candidates that were expected to be at higher redshift, but 
also contains a larger number of lower redshift objects than 
the VLT sample. These are objects from the D3 field which 
is not visible in the southern sky. The mean redshifts illus- 
trated by the dotted line are zcem = 0.63 and ~zvlt = 0.60. 
The median of the VLT sample is similar to its mean, but 
for the Gemini sample, due to the large number of points at 
z > 0.8 has a median value of 0.70. 



3.2 Photometric Data and Host Galaxy Template 
Fitting 

For this study, SNLS multicolour photometry of the super- 
novae and their host galaxies was used for two main pur- 
poses: firstly for subtraction of the host galaxy contamina- 
tion from the spectra; and secondly, for comparison of spec- 
tral features with photometric properties of the supernovae. 
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The photometric lightcurve data for the high-z SNLS 
sample above, and the low-z comparison sample described 
in Section [4] was f itted using the SiFTO lightcurve- fitter 
l|Conlev et al.l 120081 ) to ensure homogeneity in the sample. 
This provided measurements of lightcurve stretch, colour 
and date of maximum. It also permitted the interpolation of 
supernova magnitudes in g' r' i' z' -bands to the date of spec- 
troscopic observation. 

From stacked images, the host galaxy surface bright- 
ness in u* g' r' i' z' -bands at the position of the supernova 
was measured. These were used to find a best-fitting 
ga laxy template spectrum. The same template library as 
in !e08| w as used, which was drawn from the PE GASE pro- 
gramme (|Fioc fc Rocca-Volmerangelll997 , ll999T) . and a star- 
burst template (SB1) from lCalzetti et al.l ( 1994 ) was added. 



Emission-line features were interpolated over. 

A \ 2 minimisation was used to compare synthetic pho- 
tometry from the templates in the u* g'r'i' z' -bands, ftem P ,i 
to the measured values f ga i,i by minimising 



^2 _ (fgal,i 



(AX + B)f t 



(i) 



gal,i 



where a ga i.i is the error on the galaxy flux in band i. The 
template is normalised by the function AX + B to allow for 
reddening in the host galaxy. The best-fit galaxy template 
(AX + B)ftem P is denoted below as fg*(\). 



3.3 Galaxy Subtraction Method 

Here we use a gal axy subtraction method similar to the one 
de scrib ed in IE08i . rather than the correction factor method 
in iBOSl . but with o ne large difference. The objects observed 
as part of the lE08l study were all observed at parallactic an- 
gle where supernova light-loss from the slit was minimised 
whereas in this study the observations were not taken at par- 
allactic angle. In order to carry out a galaxy subtraction, we 
assume the galaxy is a flat background and the supernova is 
modelled as a point-source with a FWHM equal to the see- 
ing in the band used for object acquisition (either r' or i'), 
and scaling with wavelength as A~ ' 2 . Under this assump- 
tion, light from the supernova will be lost due to atmospheric 
differential refraction (ADR). We account for this in Section 
13X31 

The galaxy subtraction method is summarised in the 
following sections. 



3.3.1 Correcting To Compensate For Different Observing 
Conditions 

Starting with the raw spectrum, ftot"(X), which is the sum 
of supernova flux, /saT(A), and galaxy flux, fg aw (X), the aim 
is to produce a final corrected supernova spectrum, f CO rr(X). 

The SNLS observations were not spectrophotometric, 
therefore before the galaxy template (which is normalised 
to SNLS host photometry) can be subtracted, the observed 
spectrum must be placed on the same flux scale. This is done 
by multiplying by a correction factor A sens calculated by 
comparing photometry of the supernova and host from the 
imaging, correctly scaled to match observing conditions and 
observational set-up, to the synthetic r'- and i'-band flux in 



the observed spectrum. A 3en3 is the mean of the corrections 
needed to match the r'- and i'-band flux. 

ftot(X) — ftot"(X) x A sens 

We denote the photometrically-normalised SN and host 
galaxy spectra as 

fsN(X) = fsN°(X) X A sens 

3.3.2 Template Subtraction 

The fractional contamination (FC) of galaxy light in the 
total spectrum is calculated for both the r'- and i' bands 
using the galaxy and SN photometry and accounting for the 
effects of the observing conditions. Any object with frac- 
tional contamination of galaxy light in the total spectrum 
FC ^ 0.1 in either band has the best-fitting galaxy template 
subtracted-off. The FC is calculated in two bands instead of 
one because as it is measured in the observed frame, there 
will be a redshift-dependence on the value of FC as the filters 
will fall at different points in the rest-frame of the object. In 
situations with z < 0.7, the i'-band filter falls in the redder 
part of the spectrum where there is less supernova flux rel- 
ative to galaxy; and at z > 0.7, the reverse is true. The FC 
will also be phase-dependent as the supernova spectrum be- 
comes redder with time. The value of 0.1 was chosen because 
at this level of contamination, the error on the EW measure- 
ments will be «10%, which is comparable to the size of the 
error introduced by misidentifying the galaxy host type (see 
Section[33J- Also, below this level of contamination, the risk 
of subtracting too much galaxy template is high. 

Before subtraction can occur, the best-fit galaxy tem- 
plate, fg(X), must be scaled to account for differences be- 
tween the photometry apertures and the spectroscopy aper- 
tures. The galaxy photometry is in units per square arc- 
second (or effectively in an aperture of size 1" x 1"). The 
aperture size of the spectroscopy is given by the slit width, 
which is always 0.75" for CMOS observations, but varies 
for the VLT spectra, multiplied by the size of the aperture 
window used during extraction, and is denoted by A ap in 
units of square arcseconds. Under the assumption of a flat 
host galaxy background, the flux of the host transmitted 
through the slit, f g (X), is then 

fgW = /f(A)xA aj) 

The galaxy template is then interpolated to the same wave- 
length scale as the observed spectrum and subtracted off 
leaving the intermediate step 

/corr (A) = ftot (A) — /* (A) 

= fsNW + fgW-fgW- 

In an ideal situation, where the galaxy template matches 
the observed galaxy perfectly, /| = f„ leaving a spectrum 
containing only supernova flux. 



3. 3. 3 Correcting for ADR 

The spectrum now contained only supernova light, but the 
colours were not correct due to the lo ss of light at the w ave- 
lengths not used for acquisition. See iFilippenkol (1 19821 ) for 
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Figure 2. Example of host galaxy subtraction from a high-z 
SNLS supernova spectrum. The raw data (for object 07D3ea at 
z=0.471) is shown in light grey, with a smoothed version over the 
top in black. The SN spectrum after subtraction is shown as the 
smoothed light grey line underneath. 07D3ea was best-fit by an 
star-burst galaxy template and has FC = 0.19 in the i'-band. 



the equations governing light loss as a function of wave- 
length, hour-angle, seeing, airmass etc.. A correction func- 
tion, cadr{X), for the supernova spectrum was found by 
fitting a second-order polynomial to the inverse of the frac- 
tion of SN light which is transmitted through the slit at 
the effective wavelength of each band, i.e. multiplication by 
cadr(X) should have corrected the spectrum so it repre- 
sented the total supernova flu x. The fitting was calculated 
using MPFITFUN, part of MPFIT (|Markwardtll200St ). The final, 
subtracted, ADR-corrected spectrum, f CO rr(X), was given 
by: 



f int2 {\\ 

J corr V'v 



Note that this correction was not applied to the VLT spec- 
tra. 



3.3.4 Extinction Correction 

The final step is to correct for extinction in both the host 
galaxy and the Milky Way . In the rest-frame, the SALT2 
colour law (|Guv et al.ll2007l ). a function of wavelength, red- 
shift and SN colour, is applied to each object to account for 
the host extinction. This correction is small and has almost 
no effect on the EW measurements . In th e observed frame, 
the Cardelli, Clayton and Mathis l|l989l ) extinction law is 
applied to account for Milky Way extinction using the dust 
maps of ISchleeel etHI (|l998l V The SNLS fields are selected 
to be in regions with low galactic extinction so the Milky 
Way correction is not large either. These corrections can be 
summarised as c ext (X, z, c) where c is the supernova colour. 



J corr (A) 



(/siv(A) + /g(A) ~ fg{X))cADR{\)c ex t{\z,c). 



The final spectrum contains only supernova light and can 
have the equivalent widths of the features measured. 

The same scaling treatment is applied to the error 
spectrum, but no subtractions are done because it is as- 
sumed to contain only sky noise. As an example, the galaxy- 
subtracted spectrum for 07D3ea is illustrated in Figure [2] 



Table 2. The spectral ranges defined in Folatclli (2004) compared 
to the new ranges defined here. 



Feature 


Folatelli Blue 


New Blue 


Folatelli Red 


New Red 




Limit (A) 


Limit (A) 


Limit (A) 


Limit (A) 


Call 


3500 - 3800 


3500 - 3700 


3900 - 4100 


3900 - 4100 


Sin 


3900 - 4000 


3900 - 4000 


4000 - 4150 


4050 - 4125 


Mgn 


3900 - 4150 


4050 - 4200 


4450 - 4700 


4450 - 4700 



3.4 Equivalent Width Measurements 

The equivalent widths were measured in the same way as in 
iBOSl . fitting a pseudo-continuum to the feature and calcu- 
lating the EW between these limits. The pseudo-continuum 
was defined by the highest point in a region either side of 
the feature. In general, an EW and its associated statistical 
error are given by: 



EW = 



2 

°~EW 




(2) 
(3) 



where f\ is the flux of the spectrum, / c is the flux of the 
pseudo-continuum, 07 is the uncertainty in the flux from 
the error spectrum and a Ci is the uncertainty in the pseudo- 
continuum. However, as error spectra were available for the 
SNLS observations, it was possible to weight each point in 
the calculation and account for uncertainties in the flux due 
to sky noise. Equation [2] was modified to: 



EWst- 



E 



1 - 



/c(Ai) 



AXiWi 



N 



(4) 



Wi is the weighting of each point given by Wi = 
. The statistical error was still calculated as given 



where 

in Equation [3] 

The same EW code for measuring the features as in 
lB08l was used with changes to the method used to select the 
limits of the pseudo-continuum to make the process more 
automatic. The ranges of wavelength for the blue and red 
bounds of the SNe la featu res were first defined in |F04| and 
were the ones used in IBOSl . In this study, we used the Ca 11 
H&K, Sill and Mgn features and the boundaries for these 
features were redefined according to observations of low-z 
and high-z objects using trial and error to find the optimum 
values in order to make the process more automated. The 
new ranges are given in Table [2] 

As the definitions of the SNe la spectral features in 
|F04| overlap, the code could over-estimate the Can feature 
or the Mgn feature by including the Sin feature. The code 
now automatically adjusts the ranges to compensate for this 
depending on the supernova "colour" . As a crude way of cal- 
culating this from the spectrum, the ratio of integrated flux 
in the region 50-100 A above and below the minimum of the 
Si n feature was measured and used to move the boundaries 
of the features so the over-estimation no longer occurs. For 
a redder supernova, the red limit of the Call feature was 
decreased and for a blue supernova, the blue limit of the 
Mgn feature was increased. 

In the event that the automation still failed, the bounds 
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Figure 3. A histogram plotting the minimum x 2 values found for 
each of the galaxy templates for the photometry of 06Dlbz. The 
best fit is an SO galaxy, but the elliptical and irregular templates 
cannot be ruled out. The lower line marks the minimum x 2 an d 
the higher line marks the 3<r region. There is no x 2 value for the 
Sa template as this fit does not converge. 



Another change from lBOSI is that measurements of EW 
features where a GMOS chip gap has fallen in the boundary 
regions of the feature are no longer removed from the sample. 
Instead, a linear interpolation is made across the gap. The 
presence of the chip gap adds to the error on the EW mea- 
surement as there is no data in that region. To compensate 
for this, the statistical error is increased by a multiplicative 
correction factor A gap to take into account the number of 
pixels which form the feature where there is essentially no 
information. 



A a 



Ntot 



N tot - N„ 



where Ntot is the total number of pixels in the feature and 
Ngap is the number of pixels from the chip gap included in 
that feature. 



4 LOW-Z COMPARISON SAMPLE 



for the pseudo-continuum had to be selected manually. In 
this case an error to account for the uncertainty in the place- 
ment of the feature bounds was calculated by varying the 
boundaries within ±10 A at random and evaluating the EW 
fifty times. The resulting EWs were fitted with a Gaussian 
and the pseudo-continuum contribution to the EW error was 
given by the standard deviation of the Gaussian. This was 
added in quadrature to the statistical error. 

The EW measurements for the 149 objects in the final 
sample are in Tables [A"3l I A4I and I A5I At this point, five ob- 
jects were removed from the sample as the galaxy subtrac- 
tion method was believed to have failed for these objects. 
This was reflected in either negative values of EWs or values 
of EW Mg II close to 0. These objects tended to occur in the 
very centre or on the outskirts of the host where the assump- 
tion that the host contribution does not vary within the size 
of the spectroscopy aperture. The objects are marked with 
an * in the results tables. This left a total of 144 objects: 78 
objects from the Gemini telescope and 66 from the VLT. 



3.5 Evaluation of Additional Sources of Error 

Using the multicolour host galaxy subtraction method there 
are several extra terms in the error calculation as well as 
the statistical error. Firstly, there is an uncertainty in the 
selection of the galaxy template. Any galaxy template within 
a 3cr range of the best-fit galaxy is identified to rule out 
any galaxy template which is definitely not the host. This 
is illustrated in Figure [3] for the example of 06Dlbz above. 
The best-fitting template is that of an SO galaxy, but the 
elliptical and irregular galaxy spectrum lie within the 3<r 
range. 

Any galaxy template fit falling within this 3<j range 
are appropriately scaled, subtracted from the raw spectrum 
and the EW measured. The scatter in EW measurements is 
then used to produce an error due to misidentification of the 
galaxy type. It is typically in the order of a few A (i.e. 10% 
EW) and is added in quadrature to the la statistical error 
from the EW measurement. 



4.1 Spectroscopy Sample 

The equivalent widths of the Can H&K, Sin and Mgn fea- 
tures change with the phase of the supernova so it is impor- 
tant to find a sample with good phase coverage. In order to 
select a sample, the SUSPECT archival was used to search 
for published spectra with phases of ±11 days from B-band 
maximum. This phase range was chosen to match the range 
of rest-frame phases in the high-z spectra. Before —11 days 
spectra are rare because of the time needed to discover the 
objects and arrange spectroscopy scheduling. After 11 days, 
there are not many high-z spectra as the supernovae have 
faded. The selection criteria applied to the spectra were: 
that there should be at least two spectra within the ±11 day 
range as we want to study the evolution of spectra features 
with time; the spectra should include all three of the Call 
H&K, Sill and Mgn features so as to be certain the super- 
novae are normal in the whole region being investigated; all 
the objects should be classified as core-normal supernovae. 

There are a number of problems with applying these cri- 
teria to the current low-z sample. Firstly, a lot of published 
spectra do not have the wavelength coverage to include the 
Ca II H&K feature as this is in the near-UV which is difficult 
to observe from the ground. Secondly, the well-studied ob- 
jects tended to be the spectroscopically or photometrically 
peculiar examples which are not useful for this study. 

Applying these criteria to the SUSPECT archive pro- 
duced 104 spectra from 14 objects. This wa s then combined 
with 4 9 spectra from 11 objects published in lMatheson et all 
(2008); however, the latter objects are excluded from the 
analysis of the Call H&K feature. This is because the spec- 
trograph used has a blue limit of 3700 A which would mean 
the feature would only be completely visible in spectra of 
objects with a narrow Call feature at certain redshifts and 
so bias the sample. Table details the 23 low-z objects used 
in this study. 



4 http : //bruf ord.nhn. ou. edu/" suspect/ 
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Table 3. Low-z sample used to make the low-z trends. Objects 
with published photometry are marked with a t symbol. Objects 
in galaxies with known distance moduli or within the Hubble 
flow {z ^ 0.2) are marked with a t. Five objects have been used 
in the low-z cosmology sample in lAstier et al. (2006) and these 
are marked with an *. 



Supernova 



Phases of Spec- 
tra Used 1 



Reference 



SN1989B tt 
SN1990N tt 

SN1994Dtt 

SN 1996X t 
SN1998aq tt 
SN1998bu tt 
SN 1999ee t 

SN 2000E t 
SN 2002bo t 

SN 2002er t 



SN 2003du t 



SN 2004S t 
SN 2004dt 

SN 2004eo t 
SN 1997do t 
SN 1998V t* 
SN 1998dh t 
SN1998ec tt 
SN1998eg tt* 
SN1999cc tt* 
SN 1999ej 
SN 1999gd t 
SN1999gp tt* 
SN2000fa tt* 
SN 2001V t 



Barbon et ah (1990) 
Lcibundgut et al. (1991) 
Mazzali et al. (1993) 
5Patat et al. (1996) 



0,6,11 
-8 
2 

-11,-10,-9,-8 
-4,-2,2,4.10,11 

-2,0,1,7 Salvo et al. (2001) 

-9. -8.1. 2.3.4.5.6. 7 Branch et al. (2003) 
-4,-3,-2,8,9,10.1 1Jha et al. (1998) 
-11.-9.-4.-2.0.5. Hamuv et al. (2002) 
7,9 

-9,-6,-5,5 Valentini et al. (2003 

-11.-6.-5.-4.-3, Benetti et al. (2004) 
-2,-1,4 

-7,-6. Kotak et al. (2005) 
2 



-9,-8, 



-11 

-5,-4,-3 
2,4,5,6,10 
-11,-7 

-8,-6,-4,-2 
1,3,4,7,8,9,10 
1,7 

-10,-9,-7,-i 
-3,-2,-1,2,3 
-6 - 3,2,7,11 
-11,-10,-7,- 
1,2,3 

-9,-8,-7,-5 

-3,-2 

0,5 

-3,-1,1,2 

-1,3,5,9 

3,10 

-5,-2,1,3,5,7 
-10,2,3,5,10 
-11,-10,-9,- 
-6,-4,9,10 



-1,0, 

Anupama, Sahu & Jose (2005) 
-l.O JStanishev et al. (2007) 

j£risc^unas ii et_al i (2007) 
-4. Altavilla et al. (2007) 



Pastorcllo et al. 
6 , 9 , 1 M at heson_et i ^ih i 

Mathesonetah 

-3.0 Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 
Matheson et al. (2008) 

8.- 7Matheson et al. (2008) 



4.2 Photometry Sample 

The photometric properties of all 21 supernovae which have 
published photometry were refitted using the SiFTO light 
curve fitter to produce a consistently-fitted sample. The 
lightcurves were fitted for stretch, colour, date of maximum 
and peak magnitude, which were useful for studying any pos- 
sible relation between the spectroscopy features and residu- 
als from the Hubble diagram. The results are summarised in 
Table lA"6l which also shows any independent distance mod- 
ulus measurements available for the host. 



4.3 Low-z Equivalent Width Measurements and 
Errors 

The EWs were evaluated using Equation [2] and all measure- 
ments were made in the object's rest-frame. The values are 
given in Table IA7I Error spectra were not available for the 
low-z spectra as they were for the high-z SNLS data which 



is why Equation [2] was used. This meant that no statisti- 
cal e rror on the EW could be measured. The assumption in 
|B08| was that the main source of error was due to incorrectly 
placing the pseudo-continuum bounds around the features 
and this error was quantified in the same way as described 
for the high-z data in Section [3.41 however, the low-z spec- 

-tra are of high signal-to-noise and, as the equivalent width 
measurement code had been made much more automated, 

_these errors were negligible. When one measurement of the 
EW was required per supernova, the mean over all epochs 
was measured and the error was assigned based on the RMS 
scatter about the mean. This was a good estimation for the 
features being investigated here as the spectral features ap- 
pear to be approximately constant with phase within the 
phase range investigated (see Figure \%$. 

4.4 Low-Redshift Trends 

We calculate low-redshift trends as a measure of the varia- 
tion of spectral features with phase relative to the date of 
maximum. The trend is calculated by measuring a resistant 
median each day ±10 from B-band maximum. The resis- 
tant median removes outliers which are more than 3<r from 
the median. The weights applied to each measurement is a 
Gaussian centered on that date with a FWHM of 4 days. 
This time-frame was chosen as smaller time intervals result 
in the trend becoming affected by the scatter of the EW 
and any poor phase sampling. A weighted RMS was calcu- 
lated to produce the lcr scatter for the sample. The resulting 
trends are presented in Figure 3] 

The trends shown in Figure U show that over the range 
of phases studied, all three features are approximately con- 
stant. The Call H&K feature decreases in strength slightly 
post-maximum, but only by 20 A. Both the Call H&K and 
Mgn features show one supernova with a significantly re- 
duced EW measurement. In the case of Call H&K, it is 
SN 1998aq and for the Mg II feature is it SN 1994D. Both ob- 
jects appear to have quite a blue colour (see Table 1X6)) . Due 
to the small size of the sample, especially for the Call H&K 
feature, it is not possible to tell whether these are anomalous 
results, or important sub-classes of objects which have not 
yet been identified. For example, if one object in the low-z 
sample has a low Ca II H&K equivalent width then a number 
of objects like this should these should appear in the high-z 
sample. This may explain some of the low EW values seen 
in Figure [5] 

Thes e tre nds differ in a number of ways from those pre- 
sented in |B08| . Firstly, the metho d of calculating the trends 
is different. The trends in |B08| a re created by measu ring 
EWs from the template spectra of iNugent et al.l (|2002l ). An 
error region was determined by adding an error to the trend 
and comparing to low-z measurements until a Xr — 1 was 
reached. This method assumes that the low-z measurements 
are distributed symmetrically about the trend. 

In this paper, the trends are generated directly from 
the measurements themselves and the la scatter was deter- 
mined directly from the EW measurements. There are also 
a numb er of objects included in this study which were not 
in iBOSi and vice versa due to the selection criteria or date 
of publication. As the sample size is small, this has had a 
consider able e ffect on the trends. 

The IBOSi Si II trend is slightly lower that the one seen 
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low-redshift trends seen in Figure [4] with the high-z points 
measured from the Gemini and VLT samples overplotted. 

The high-z measurements for the Si II feature appear to 
be smaller than the low-z trend. In order to look for system- 
atic shifts, we make trends of the high-z data. The trends 
are calculated in a similar way to the low-z ones, except that 
the points are weighted according to the size of the error on 
the measurement, rather than applying a Gaussian weight 
based on phase. A rolling mean is applied with a bin-size of 
±3 days. The mean for each phase bin is constrained so that 
there must be more than 10 equivalent width measurements 
within ±3 days in order for the value to stand. Consequen- 
tially, the trends are not evaluated at some early and later 
phases. 

Figure [5] shows that for the Ca II H&K and Mg II fea- 
tures, the low- and high-z trends are approximately consis- 
tent. The Sill feature shows a high-z trend which is consis- 
tently below the low-z one, as expected from the low-EW 
points seen in Figure [5] This will be discussed later. This is 
contrary to the observations of Bo3 which saw no difference 
between the high- and low-z populations in the Sill feature. 
The difference between the high- and low-z populations that 
may have been observed for the Mgn feature in B08 is not 
seen in this analysis which uses a larger high-z sample and 
a different method for generating low-z trends. 



Figure 4. Low-z Trends for Can H&K, Sill and Mgn. The indi- 
vidual supernova measurements are marked by the different sym- 
bols. The central thick dotted line is the trend and the grey shaded 
region represents the la error distribution. The legend applies to 
all 3 panels. 



here. This is primarily caused by the fact that new spectra 
intr oduc ed into the sample all had higher EW values than 
the |B08| mean. The increase in EW to maximum followed 
by a decline is not seen here to such a large degree as in the 
IB08i trends. 

The rapid rise in the Mg II EW seen after +3 days in BOS 
is not seen in this sample where the trend is flat from —10 
to +10 days. This is likely due to the fact that the sample 
is different; however, SN 1989B, which is in both samples, 
does show a larger increase in the Mgn EW between the 
measurements at +6 and +11 days. 

Due to the lack of available spectra of SN 1991T-like 
and SN 1991bg-like objects, it was not possible to measure 
trends for these objects in a similar way to the core-normal 
supernovae, which is why spectroscopically peculiar events 
are removed from the EW analysis. 



5 RESULTS 

5.1 Comparison between High-z Data and Low-z 
Trends 

The high-z sample is divided into the Gemini sample (ob- 
ject s pre sented in this paper, plus those published in |H05| 
and IB08T ) and VLT samples in order to identify any system- 
atic differences between the samples: none are seen in this 
analysis. Figure [5] shows the comparison for the three fea- 
tures in this analysis for both samples. The figures show the 



5.2 Correlations Between Spectroscopic and 
Photometric Properties 

One of the reasons why there might be differences between 
the low- and high-z samples is because the spectral prop- 
erties will be sensitive to the explosion physics of the su- 
pernovae. There is increasing e vidence that t h ere a re two 
populations of supernova (e.g. ISullivan et al.l (|2006')) and 
that the prompt supernovae tend to be brighter and bluer. 
This will produce a drift in the mean photometric proper- 
ties of SNe la as a function of redshift if the mix of prompt 
and d elayed supernovae changes with redshift l|Howell et al.l 
2007). Mb , the absolute B-band magnitude, is known to cor- 
relate with lightcurve width in the sense that more luminous 
SNe have higher stretch. We also know from |B08| that EW 
Si II is correlated with M b , so brighter supernovae appear to 
contain less Sill. These facts taken together imply that the 
mean spectral properties between the high- and low-z sam- 
ples will be different in the sense that the high-z sample, 
where we expect a greater proportion of the prompt-type 
objects, should have a lower EW Sin. This is indeed the ef- 
fect seen in Figure [5] In the following we quantify the size of 
this effect and determine whether it explains the observed 
differences seen in Figure [5] First we quantify the correla- 
tions between spectral and photometric properties, stretch 
and colour, using a Spearman's Rank correlation method. 

We combine the Gemini and VLT samples into one high- 
z sample as the correlations are likely to be small effects. The 
high-z spectra in this study are also noisy and so there are 
large errors on the EW measurements. As the correlation 
coefficient calculation does not deweight objects based on 
the size of their error bars, we have grouped all the high- 
z objects together and ranked the objects according to the 
size of their error and used only the lowest 75% for this 
part of this analysis. The cuts for Call H&K, Sill and Mgn 
correspond to error values of 18.1 A, 9.7 A and 19.3 A re- 
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Figure 5. Each figure shows the low-z trend of the EW of either Ca II H&K, Si II or Mg II as a solid black line with the la region in grey. 
The high-z sample is shown as the points plotted over the top as squares and circles (Gemini) and triangles (VLT) 



spectively. The correlations are evaluated for all the high-z 
objects with errors less than the 75-percentile cut. The cor- 
relation coefficients are presented in Table [3] 



The correlation analysis shows that the strongest corre- 
lations exist between stretch and EW Sin, and this is seen 
in all samples at the 2.5-3.5rj level. There are also weak 
correlations between the EW Sill and colour for the com- 
bined sample, and for the low-z sample between EW Mgn 
and colour. There also appears to be a slight correlation 
between the EW Call H&K and colour. There appears to 
be no correlation between EW Call H&K and stretch. We 
now examine the two most promising relations: EW Sill and 
stretch; and EW Mgn and colour. 



5.2.1 EW Sill and Stretch 

The stretch-EW Sill measurements are plotted in Figure [7] 
for the low-z (diamond) and high-z (square and triangle) 
points. Here the high-z Gemini and VLT samples are plot- 
ted separately in order to identify any possibly systematic 
offsets between the samples: none are seen. The best-fit was 
found using LINFITEX, part of MPFIT, which takes account 
of errors on both the variables. The low- and high-z rela- 
tionships shown by the black lines are are: 

s = (1.22 ± 0.03) + (-0.015 ± 0.002) x EWsui low<£) 
s = (1.21 ±0.01) + (-0.017 ± 0.001) x EWsm higH® 

The reduced \ 2 f° r * ne fits are Xr = 2.2 (low-z) and \r = 
2.7 (high-z). The fits are consistent to within la, although 
the errors on the gradient and intercept are correlated. We 
would expect the fits to be the same if the high- and low-z 
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Figure 6. Each figure show the low-z trend of the EW of either Can H&K, Sin or Mgll as a solid black line with the lc region in grey. 
The high-z trends are overplotted in blue (Gemini) and green (VLT) with the shaded region representing the lcr RMS scatter 



populations behave in the same way; however this does not 
mean that they must have the same distribution at high and 
low redshift. 

It has already been established in iBOSl that there is a 
relationship between EW Si II and Mb , the absolute magni- 
tudes of the individual supernovae: the brighter the super- 
nova, the smaller EW Sin. This has been extended here to 
use lightcurve stretch instead of magnitude, as is shown in 
Figure as Mb is known to correlate with stretch. This 
relationship can potentially explain the difference seen be- 
tween the low-z and high-z objects in Figures[5]and[6]above. 
Figure [8] shows that the stretch distribution is shifted to 
higher stretch for the higher redshift samples compared to 
the low-z one. The relationships seen in Figure [Jj show that 
more higher stretch objects would mean more low-EW Sill 
measurements. 

To test whether the difference in mean stretch between 



the high- and low-z samples can be attributed as the cause 
of all the systematic shift in EW Si II as seen in Figure [5] 
a predicted high-z mean EW was calculated based on the 
mean high-z stretch and the low-z relation from Equation 
(5). This predicted high-z mean EW Sill can then be com- 
pared to the measured value. 

From Figure[8] the mean stretches for the high- and low- 
z samples are si ow _ z = 0.96, SQ em = 1.00 and syLT = 1.00, 
with standard deviations o"i ow _ z = 0.10, o"Gem = 0.11 and 
°"VLT ~ 0.10. Using the high-redshift standard deviation to 
produce upper and lower bounds and ignoring the errors on 
the fit in Equation (JSJ), the predicted mean high-z EW Sill 
values are 14.67 A with a standard deviation of 7.34 A (Gem- 
ini), and 14.67 A with a standard deviation of 6.67 A (VLT). 
Taking the weighted mean and standard deviations of the 
observed high-z trends in Figure[6]yields observed mean EW 
Si II of 11.23 A with a standard deviation of 4.58 A (Gemini) 
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Table 4. Correlation coefficients between the measured EWs and 
the light curve stretch and colour for the low-z, high-z and com- 
bined samples. The number of standard deviations from a null 
hypothesis is given in the final column. 



Feature 


Variable 


Sample 


P 


a from null 


Ca IT H&^K 


Stretch 


T 


0.21 


-0.66 








0.10 


-1.0 






All 


0.09 


-0.91 




Colour 


Low-z 


0.44 


-1.38 






High-z 


0.12 


-1.20 






All 


0.21 


-2.15 


Sin 


Stretch 


Low-z 


-0.76 


3.58 






High-z 


-0.24 


2.50 






All 


-0.33 


3.78 




Colour 


Low-z 


0.19 


-0.89 






High-z 


0.08 


-0.80 






All 


0.19 


-2.12 


Mgn 


Stretch 


Low-z 


-0.31 


1.47 






High-z 


0.09 


-0.93 






All 


0.03 


-0.38 




Colour 


Low-z 


0.47 


-2.22 






High-z 


0.09 


-0.85 






All 


0.17 


-1.87 



and 11.05 with a standard deviation of 6.64 A (VLT). Al- 
though the errors on the quantities are large, the predicted 
and observed high-z EW Sill are in good agreement with 
each other. 

Thus it appears that the difference in observed EW Si n 
between the low- and high-z samples can be explained by the 
difference in the stretch distributions of the two samples, 
i.e. the differences are due to changing demographics rather 
than a new, unexplained form of evolution. 
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5.2.2 EW Mg II and Colour 

The relationships between colour and EW Mgn in the low- 
and high-z data are plotted in Figure [5] Figure highlights 
the differences in properties between the high- and low-z 
samples. The range of EW Mgn values in the high-z data is 
much broader than for the low-z sample, but the reason is 
unclear. This could be because the low-z sample is too small 
to be representative of the SNe la population as a whole. 
The high-z sample is suffering from selection bias as a rough 
colour-cut is applied during high-z target selection and we 
do not see the range of supernova colours observed in the 
low-z sample. 

Figure [10] shows that the colour distributions for all 
three samples are very different. The low-z sample contains 
a number of redder objects which aren't seen in the high- 
z data. These are not underluminous, redder supernovae, 
but objects classed as core-normal. The Gemini distribution 
shows a strongly-peaked distribution around c = whereas 
the VLT distribution is more spread out. The mean colours 
for the samples are ci ow _ z = +0.08, CQ em ; m = —0.03 and 
E VLT = -0.04. 

It is not possible to state anything conclusive about 
any relationship between colour and EW Mgn as the data 
is very noisy at high-z and the sample at low-z is very small. 
This would be interesting to investigate further with better 
data as the relation could potentially reveal more about the 



Figure 7. Top: Stretch against EW Si II for the low-z data. Bot- 
tom: Stretch against EW Si II for the high-z Gemini (blue squares) 
and VLT (red triangles) data. In both plots, the black line rep- 
resents the best-fit straight line to the data. The fits are given in 
Equations J5J and iftj) 

supernova progenitor and be another test for an evolving 
population, but it is not possible at this time. 

5.3 Applications of SNe la Spectra in 
Cosmological Studies 

One of the aims of this study is to investigate potential uses 
for supernova spectra in the calibration of the objects for 
cosmological studies. This could be either as an alternative 
to the current method of photometric stretch and colour, or 
as a third calibrator, with the goal of reducing the scatter 
around the best-fit line on the Hubble diagram and hence 
improve the accuracy of cosmological parameter estimation. 
The large, spectroscopic dataset, with accompanying pho- 
tometry lends itself well to this type of study. 

5.3.1 Cosmology Fitter 

In order to measure a scatter in the Hubble diagram, we 
must first fit the cosmology. In this section we detail the 
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Figure 8. Histograms of the stretch distribution for the low-z 
(green), VLT (red) and Gemini (blue) samples. The dotted line 
marks the mean stretch for each sample. The low-z sample has a 
lower mean stretch than the two high-z samples. 



cosmological fitter used for this purpose. Throughout, we 
assume a flat universe with w = — 1. 

The supernova cosmology is fitted using a MINUIT- 
based, lighter- weight version of SIMPL E COSFITTEij! ] as op - 
posed to a full Bayesian analysis used in lAstier et al. (|2006t ). 
We refit the data r ather than using the cosmology found in 
lAstier et alJ l|2006l ) because the samples used in the stud- 
ies are different. A best-fit method is faster and sufficiently 
precise for this study. 

SIMPLE COSFITTER takes the supernova apparent B- 
band magnitudes at maximum, m_g; the redshift, z; the 
supernova lightcurve stretch, s; and supernova lightcurve 
colour c and performs a least-squared fit to find the best- 
fitting values of CIm, a, /3 and M by minimising the x 2 
where 



(M-51og(d L ( z ,6»)/10pc)) 2 



objects 



2 

in t 



and 



ix = tub — M + a(s — 1) — /3c 



(7) 



(8) 



9 stands for the cosmological parameters that define the fit- 
ted model, which in this case is simply Qm', c(aO i s the 
error on the observed distance modulus; and a int is the 
intrinsic dispersion in the supernova absolute magnitudes 
|Astier et al.ll2006r i. A good fit is obtained by varying the 
parameter a int manually until the Xr — lj the other param- 
eters are varied by COSFITTER itself. By assuming a value of 
Ho then the mean supernova absolute magnitude Mb can 
be calculated from the fitted value of M. 



http : //qold. astro .utoronto . ca/conley/simple_cosf itter/ 
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Figure 9. Top: Colour against EW Mgll for the low-z data. Bot- 
tom: Colour against EW Mg II for the high-z Gemini (square) and 
VLT (triangle) data. 



The code can be exploited to fit alternatives to stretch 
and colour, e.g. equivalent width. It was also adapted to fit 
for a fourth variable which has a linear contribution to the 
distance modulus. For instance, the fit could be performed 
using equivalent width as a third possible calibrator: 



/'• 



■M + a(s-l) -Pc + jEW 



(9) 



The code also enables evaluation of .M, a, f3 and 7 within a 
fixed cosmology. 



5.3.2 Cosmological Parameters Derived From 
Photometric Data 

In this section, we fit the cosmology to the data sample in 
order to have a baseline to compare against when testing for 
spectroscopic calibrators. 

As stated in Section [4.21 the low-z sample used for the 
spectroscopy study is not sufficient for cosmology-fitting. 
This is because most of the objects are not in the Hub- 
ble flow and their peculiar velocities affect their redshifts 
and hence their distances are very poorly con strained. To 
account for this, we use the low-z sample from lAstier et al.l 
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Figure 10. Histograms of the colour distribution for the low-z 
(green), VLT (red) and Gemini (blue) samples. The dotted line 
marks the mean stretch for each sample. The low-z sample has a 
redder mean than the two high-z samples. 



(2006) with photometry fitted using SiFTO, but removing 
the objects marked with * in Tabled to keep the cosmology 
and spectroscopy samples entirely separate. We also make 
stretch and colour cuts on the sample to remove any possi- 
ble outliers which may not be fully described by the width- 
brightness or colour relations. The cuts are 0.75 ^ s ^ 1.25 
and —0.25 ^ c ^ 0.25 which leaves a sample of 32 low-z 
objects, 75 Gemini objects and 66 VLT objects. This is sim- 
ilar to the cuts that would be applied for a full cosmology 
analysis. 

For the first fit no stretch or colour correction was ap- 
plied to the data. This is designated Fit A in Table [5] The 
fit shows that the intrinsic dispersion (aw) must be set to 
be quite high (0.255 mag) in order to obtain a Xr ~ L The 
value of Qm = 0.462 is higher than current measurements 
which employ stretch and colour correction (e.g. lAstier et al.l 
(2006)), but this is to be expected because at higher redshift 
the supernovae tend to be brighter and bluer which causes 
the fit to favour a higher matter fraction. 

In Fit B, stretch and colour corrections are applied 
and the resulting best-fit value of S1m is much lower. The 
value is within 2 s tandard deviations of that published in 
lAstier et all (|2006T ). The latest full analysis of SNLS data 
finds, Q,m = 0.268 ± 0.020 (Sullivan et al. (in prep.)), and 
so the value calculated in Fit B is within 3 standard devia- 
tions. The value of a is very similar to the one achieved from 
the full a nalysis, and 8 is ver y different from the value pub- 
lished in lAstier et all {2006), but in close agreement with 
the latest value from the SNLS analysis (B = 2.9). This is 
due to the fact that the treatmen t of the B paramete r in the 
fitting methods has changed since lAstier et al] |2006t ). Fit B 
shows a reduced RMS compared to Fit A which shows that 
the calibration applied using stretch and colour corrections 
are important. Applying these corrections also reduces the 



Figure 11. Hubble diagram for Fit B (solid line) which uses 
stret ch and colou r corre ctions, compared to the cosmology fit 
from lAstier et al. | ||2006|. dashed). The da ta used are the low-z 
cosmology sample from lAstier et all |200fJ) , and all the VLT and 
Gemini objects from the samples defined in Section[2] Stretch and 
colour cuts are applied. In both plots, the black points are calcu- 
lated using the distance moduli from the apparent magnitude of 
the supernova fitted by SiFTO and the cosmological parameters 
from Fit B, correcting for stretch and colour. The lines represent 
the distance moduli calculated from Equation 1(5} . The bottom 
plot shows the residuals from Fit B (see Section 15.3.31 for defini- 
tion) compared to the Astier cosmology. 



intrinsic scatter required for a good fit. It is a result similar 
to this which is the goal of this part of the study. 



5.3.3 Hubble Residual Correlations 

Now the cosmology has been fitted to the data, residuals 
from the fit can be plotted against spectral features and 
possible alternative calibrators found. The Hubble residuals, 
R, as shown in Figure ITT1 are defined as 



R = m B — m CO smo(z, fit 



(10) 



where ms is the observed B-band apparent peak magnitude 
and mcosmo is derived from the cosmological parameters. 
The assumed cosmology is that measured using Fit A i.e. 
fijvf = 0.462 in a flat universe. The correlations with the 
photometric quantities stretch and colour are also presented 
as a guide as, naturally, strong correlations are expected for 
these variables. As in Section 15.21 we use only the high-z 
EW measurements from the least noisy 75% of the sample. 

For the low-z sample, the method of calculating the ob- 
served distance modulus is slightly different as not all the 
objects are in the Hubble flow. For five very nearby objects 
(see Table |A"6]) which are not in the Hubble flow we use dis- 
tance moduli measured from their host g alaxies either using 



Cepheid variable stars (ISaha et alJ 



brightness fluctuations (|Aihar et al 



2006) or galaxy surface 



20011 ). These values mi- 
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Table 5. Cosmology fits table. The fit from lAstier et al 1 l|2006l l is included for reference. Fits C, D and E were made 
without the low-z sample and fixing the cosmology to the values from Fit B. Fit F was made without the use of 
the low-z sample, but the cosmology was not fixed. Parameters which were fixed are marked with a *. The quoted 
values of <r; n t are those used to obtain the best fits. 



Fit 


Extra Calibrator 




RMS 


f2 A f 


a 


P 


7 


M 






(mag) 


(mag) 










(mag) 


Astier 








0.263±0.053 


1.52±0.14 


1.57±0.15 




23.85±0.03 


A 




0.255 


0.260 


0.462±0.065 








24.091±0.045 


13 




0.12 


0.181 


0.197±0.031 


1.28±0.14 


3.21±0.19 




23.999±0.025 


C 


Can H&K 


0.10 


0.176 


0.197* 


1.28* 


3.21* 


O.OOliO.001 


24.050±0.028 


D 


Sill 


0.10 


0.183 


0.197* 


1.28* 


3.21* 


0.006±0.002 


24.084±0.030 


E 


Mgn 


0.10 


0.177 


0.197* 


1.28* 


3.21* 


0.000±0.001 


24.044±0.039 


F 


Sin 


0.10 


0.193 


0.214±0.053 


1.50±0.18 


3.49±0.24 


0.007±0.002 


24.123±0.054 



Table 6. Correlations between photometric and spectral proper- 
ties and the Hubble residuals. There are not enough low-z objects 
with EW Call H&K measurements to measure a correlation. 



Variable 


Sample 


P 


(7 from null 


Stretch 


Low-z 


-0.45 


1.62 




Gemini 


-0.45 


4.04 




VLT 


-0.51 


4.10 


Colour 


Low-z 


0.79 


-2.73 




Gemini 


0.51 


-4.48 




VLT 


0.58 


-4.67 


EW Can H&K 


Low-z 








Gemini 


0.09 


-0.60 




VLT 


-0.13 


0.97 


EW Si II 


Low-z 


0.72 


-2.59 




Gemini 


0.13 


-0.94 




VLT 


0.10 


-0.71 


EW Mgn 


Low-z 


0.78 


-2.80 




Gemini 


0.05 


-0.38 




VLT 


-0.29 


1.87 



nus the value of Mb determined from the best-fit value of 
M from Fit A and assuming Ho = 70kms Mpc -1 , are 
used as m C osmo- For 5 other objects, they are situated in 
host galaxies with z ^ 0.02 and can be assumed to be in the 
Hubble flow and so the same method as the high-z sample 
to calculate mcosmo is used. Plots of the low-z and high-z 
data are shown in Figure [L2l The correlations are calculated 
and tabulated in the same way as Tableland presented in 
Table H 

The strong correlations between stretch and colour and 
Hubble residual in Table [6] show that the current low- and 
high-z spectral features should not be used as direct replace- 
ments for either of the photometric properties as the rela- 
tions are much weaker. The strongest correlation between 
one of the high-z spectral features and the residuals is for 
EW Mgn for the VLT sample, but this is not seen in the 
Gemini data - the correlation goes in the opposite direction. 
However, the linear fits to the data in Figure[l2]where points 
are weighted according to their errors show the best-fit lines 
all have positive gradients, implying that the negative corre- 
lation coefficient seen in the VLT data is spurious. By com- 
paring Table [5] and Figure [T2l we can see that even applying 
a cut to remove the most noisy points from the high-z data, 
the correlation coefficient does not necessarily represent the 



best way to look for a relation between two quantities as it 
is affected by outliers. 

Instead of looking to replace stretch or colour with one 
of the EW measurements, it might be more instructive to 
look for a correlation which exists between the EW and the 
Hubble residual after the application of stretch and colour 
correction. To do this, we fix the cosmology and values of a 
and j3 to Q M = 0.197, a = 1.28 and /3 = 3.21 as in Fit B, 
and use the adapted version of CQSFITTER to fit Equation 

©■ 

To see just how great an improvement can be made, 
the instrinsic dispersion was allowed to vary and Equation 
([9]) fitted allowing M, ow and 7 as the free parameters. 
The cosmology must be fixed as the objects in the low-z 
cosmology sample do not have accompanying spectra and so 
could not be used to anchor the Hubble diagram for these 
fits. The results of this process are given by Fits C, D and 
E in Tabled 

The results in Table [5] show that by adding in a third 
calibrator, the intrinsic dispersion necessary for a good fit 
can be decreased, as expected. The Can H&K and Mgn 
features show values of 7 consistent with and so they are 
ruled out as possible calibrators. However, for the Sill fea- 
ture 7 is three standard deviations from zero and produces 
a value of M consistent with the one in Fit B. The RMS is 
comparable to the one obtained in Fit B as well. 

In order to test whether this result is meaningful, the 
constraints on a, /) and Qm were relaxed so they were free 
parameters. Even without a low-z sample to anchor the Hub- 
ble diagram at low redshift, the fit using EW Sill as a third 
calibrator still gave reasonable values for the parameters, 
denoted as Fit F in Table [5] However, the RMS scatter was 
increased slightly to 0.193 mag. 

From this is is clear that Sill could have a role to play 
in SNe la calibration for cosmology. However, with the large 
errors on the high-z EW measurements, the RMS scatter 
around the best-fit cosmology do not improve and it is not 
possible to draw any firm conclusions as to the value of 7. 
Instead, in Section[5]4] we demonstrate that the relationship 
between EW Sill and stretch in the low-z sample can be 
exploited to reduce the residuals on the Hubble diagram by 
reducing the scatter on Mb measurements. 
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Figure 12. Hubble Residual - EW Plots for the low-z (left column) and the high-z (right column) objects. The best-fits to the lines are 
also shown. Solid line - low-z; dot-dashed - Gemini; dashed - VLT. 



5.4 EW Sin as a Cosmological Calibrator 

Drawing on the fact that we have seen that EW Sill and 
photometric stretch are correlated, and that there is a strong 
correlation between the Hubble residuals and the low-z EW 
Sin, we will explore the possible use of EW Sin as a calibra- 
tor for the cosmology in place of stretch for the low-z sample. 
Here the errors on the EW are smaller so it is possible to be 
more definite about the potential of the calibrator. 

We divide the objects in the low-z spectroscopy sample 
which have photometry into two independent sub-sets in or- 
der to have a sample to define the relation and a separate 
one to test its use. The first sample (i) contains the objects 
with no reliable absolute distance measures, i.e. those with- 
out Cepheid or surface brightness fluctuation distances and 
not in the Hubble flow (those marked with a ' only in Ta- 
ble^; and the second sample (ii) containing the rest (those 
marked as The two samples are illustrated in Figure 
1131 as open circles for the objects with distance measure- 



ments and closed circles for those objects without. Using 
sample (i), the relation between stretch and EW Sill is re- 
fitted to derive an expression for "spectral stretch", s S pec, 
using LINFITEX. The best-fit is given by 

sspec = (-0.021 ± 0.005)EW SiU + (1.29 ± 0.07) (11) 

where Xr = 1-55. This fit is different from the ones in Equa- 
tions (JSJ and © as fewer objects are used. The remaining 
10 low-z objects were corrected for stretch and also spectral 
stretch according to the relation 

M B = m B — H + a(s - 1) 

where s can either be the normal photometric stretch or 
spectral stretch as fitted in Equation (|lip . We then plot 
Mb against /i for illustration purposes, and measure the 
weighted-mean, Mb, and the RMS scatter about the mean. 
This is shown in Figure [TU The values of Mb and the RMS 
are shown on each plot and summarised in Table [7] For the 
purposes of this exercise we take the values of Qm and a 
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Figure 14. Comparing spectral stretch correction to photometric stretch correction with and without colour correction. The value of Mb 
given on the plot is the weighted mean of the points after the stretch correction. Without colour correction (left column) the RMS before 
any stretch correction is 0.45 mag and the mean absolute magnitude Mb = —18.94. With colour correction but no stretch correction 
(right column), Mb = -19.19 mag and the RMS is 0.24mag 




Figure 13. The two sub-samples of the low-z data are shown. 
Those with closed circles do not have known distance moduli 
and are used to define the spectral-stretch relations (black line). 
Those with open circles will be used to test the relation and are 
not included in the fit. 

from Fit B. The right-hand column of Figure [141 shows the 
effect of also adding a correction for supernova colour. 

Table shows that the lowest RMS is obtained when 
using photometric stretch and colour together, but when no 
colour correction is made, spectral stretch correction pro- 
duces a lower RMS than photometric stretch. 

It is not possible to apply this correction to the high-z 
supernovae and see any improvement in the data due to the 
large errors on EW Sill, but in the future with a sample 



Table 7. A summary of the vaules of Mb and the RMS scatter 
about the mean for the fits in Figure H4l 



Corrections Applied 


M B 


RMS 




(mag) 


(mag) 


None 


-18.94 


0.45 


Spectral stretch 


-19.11 


0.38 


Photometric stretch 


-19.00 


0.44 


Colour 


-19.19 


0.24 


Colour + Spectral stretch 


-19.29 


0.33 


Colour + Photometric stretch 


-19.25 


0.20 



of well-measured equivalent widths it is possible that spec- 
tral stretch could replace photometric stretch. This would 
mean that once the supernova was identified fewer photom- 
etry points would be needed to fit the light curve, only the 
B-band magnitude at maximum and the redshift would be 
needed: supernova colour could be measured from synthetic 
photometry of a spectrum obtained around maximum light. 
However, the stretch-EW Sin relation needs better calibra- 
tion first using, for example, a large sample of low-z super- 
novae either with known host galaxy distances or within the 
Hubble flow. As seen from the 7 parameter analysis in the 
noisy high-z data, as shown in Table[5j if corrections for pho- 
tometric stretch, colour and spectral stretch/EW Sin can 
all be combined then there is the potential that calibration 
could be improved. 
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6 DISCUSSION 

6.1 Tests for Evolution 

The comparison between the high-z and low-z trends seen in 
Figure [6] shows that for the Call H&K and Mgn features, 
the high-z sample is consistent with the low-z trends, but 
the scatter on the data is quite large. For the Sin feature, 
a large number of low equivalent width measurements are 
seen in the high-z data and the high-z trend is consistently 
below the low-z one. This can be explained as a change in 
the demographics. In Section 15.2.11 we see that light curve 
stretch is correlated with EW Si II. At higher redshift, due to 
the "A+B" model we expect more of the prompt component 
which are brighter and bluer as more 56 Ni is thought to be 
produced during the explosion. This would produce less IME 
which is reflected in the lower EW Sill. 

To the level of current precision, the variation in mean 
EW Si n between high and low redshift can be understood in 
terms of the combination of a change of demographics within 
a two-component SN model and an observed correlation be- 
tween photometric stretch and EW Si II. Therefore, we see no 
evidence for an additional bias that is not already accounted 
for by the standard photometric stretch and colour correc- 
tions. However, in future experiments which may go to even 
higher redshift to look for w = w(z) (models where dark 
energy is not a cosmological constant) will need a higher 
level of precision. Future cosmological studies observing to 
z = 1.7 will need a precision of 2% in the standardisation of 
supernova magnitudes if they are to keep systematic errors 
smaller than statistical errors, and it not yet clear whether 
the current method s of calibration will be suitable for this 
|Howell et al.ll2007h . We also know that the width-brightness 
relation may not be valid for supernovae with large stretch 
such as SN 2001ay-like objects. 

It will only be possible to make a strong statement re- 
garding the evolution of SNe la when there is a large, un- 
biased sample of low-z spectra. In Figure [5] we have only 3 
low-z objects in the s = 1 bin, compared to 14 in the VLT 
sample and 17 in the Gemini sample. Future surveys such as 
Palo mar Transient Fac tory (|Law et al.l 2009) and SkyMap- 
per (|Keller et al.ll2007T ) will largely address this problem. 



stretch and not as well-understood. At high redshift, colours 
may be affected by the environment surrounding the super- 
nova simply due to the fact that it is not possible to sep- 
arate the supernova light from the surrounding i nterstellar 
and circumstellar medium fsee lConlev et af1l2007i . for a dis- 
cussion of the effect of dust and the colour parameter on 
cosmology). In addition those effects, from studies in the lo- 
cal universe, there is kno wn to be an in trinsic difference in 
supernova colours as well (|Phillips!ll993T >. In Table © we see 
that supernova colour has a stronger correlation with Hub- 
ble residual than stretch, implying that it is more important 
for calibration purposes. 

Another application for spectral stretch is that it may 
allow more supernovae to be added to the low-z sample for 
cosmology. Currently, most SNe la have either good pho- 
tometric coverage or good spectroscopic coverage, but not 
both. For instance, in the spectroscopic sample used f or thi s 
study, only 5 objects overlap with the lAstier et al.l (|2006l ) 
cosmology sample. By taking the objects that have good 
spectroscopic coverage, instead of using photometrically- 
derived properties to calibrate them, spectral stretch could 
be used. This would also require a B-band magnitude at 
maximum, or interpolated from other data, and a known 
distance modulus or redshift. There is not yet a spectral re- 
placement for colour, but the EW Mg il-colour relation seen 
in the low-z sample could be exploited for this purpose when 
the relationship is better quantified. 

For future studies out to z = 1.7, iHowell et all (|2007n 
have shown that supernova calibration needs to be better 
than 2% in magnitude in order to keep systematic errors be- 
low the level of statistical errors. By introducing a third cali- 
brator based on EW Sin into the cosmology in Section [5. 3. 3 1 
we find a non-zero value for the constant 7 whilst achieving a 
similar RMS scatter around the best-fit cosmology as when 
using standard photometric stretch and colour corrections. 
The estimate of Qm is also closer to the current value from 
the full Bayesian analysis of SNLS data. In the future with 
higher signal-to- noise spectra, the EW Sin may provide a 
further method of cosmological calibration. At the moment, 
the best fit is not approaching the 2% level. This can be 
done with more observations with high signal-to-noise spec- 
tra and good photometric follow-up of discoveries. 



6.2 Implications for Cosmology 

The origin of the width-luminosity relation is well- 
understood to first order: a variation in the amount of 56 Ni 
alters the amount of gamma-rays emitted due to radioactive 
decay to 56 Co and 56 Fe which are absorbed and thermally 
reemitted by the expanding ejecta l|Colgate fe McKeell 19691 ; 
lArnett ll979Tl . The more 56 Ni produced, the wider the light 
curve " (lArnettll 19821 ). This also means the production of less 
IME, which is refle cted in Figure (see also Figure 2 in 
ISullivan et aLll2009h . What causes the varying amount of 
56 Ni is not understood although it may be due in part to 
the way that burning of carbon and oxygen occurs in the 
supernova, or the metallicity of the progenitor. The relation 
between stretch and EW Si 11 from Figure [7| is exploited in 
Section r5.4l for the low-z dataset to reduce the scatter of the 
objects on the Hubble plot by correcting Mb using spectral 
stretch. 

The origins of supernova colour are more complex than 



6.3 Improvements 

6.3.1 Galaxy Subtraction Method 

The methods presented in Section 13.31 detail the way host 
galaxy contamination is removed from the spectra. This is 
made somewhat more difficult for the Gemini objects which 
were not observed at parallactic angle and hence the spec- 
tra suffer from slit-losses. The objects from the VLT sample 
were observed using an ADC and so this is less of a problem. 
Extraction of spectra from the host galaxy and the super- 
novae separately was not usually possible due to the small 
angular size of the galaxy on the sky. In order to subtract off 
host galaxy light, we must make an assumption about the 
nature of the background light, and we assume that it is a 
flat background in the slit. This assumption is reasonable in 
most cases, but it can fail when the supernovae are situated 
in the very centre of the host, or on the outskirts. 

Another limitation is the number of galaxy template 
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spectra available for the subtraction. In order to make a bet- 
ter subtraction, ideally we would reobserve the host galaxies 
after the supernova light has faded with an identical set-up. 
However, there will be slight variations in the seeing and 
airmass so the correction would still not be perfect. An- 
other method would be to use the CFHT imaging directly 
to aid subtra ction. Photometry-A ssisted Spectral Extrac- 
tion (PHASE; iBaumont et alj 20081 has been developed and 
used on SNLS-VLT spectra to extract the supernova spec- 
trum with an extraction window deter mined from the pos i- 
tion of the target in the imaging data (|Balland et al.| [2009l. 
A ^-minimisation is used to evaluate the contribution of 
galaxy and supernova in each pixel to produce a spectrum 
composed of only supernova light. This method produces 
clean spectra for approximately 70% of observations, but 
when the supernova is close to the centre of the host, or 
the fractional contamination of the the spectrum is high, 
the method cannot separate the two components. These are 
two of the occasions when the multicolour host galaxy sub- 
traction method also fails so PHASE could be used as an 
alternative to the method used in this study, but not as a 
replacement. 



6.3.2 Low-z Trends 

In comparison to the large number of high-z objects included 
in this analysis, there are only 14 objects used to make the 
low-z trends for the Can H&K feature, and 25 for the Sill 
and Mgn features. The sample is biased towards supernovae 
in bright, elliptical hosts which will tend to host redder su- 
pernovae with lower stretch (the delayed component from 
the two-component model). This can be seen in Figure [8] 
The sample is also drawn from many sources where the pho- 
tometric coverage is different. To do a study of this kind in 
the future to an increased level of precision we would require 
a large homogeneous sample of low-z spectra and photome- 
try. 

The trends shown in Figure [4] reflect the diversity seen 
in the low-z objects, an advantage of using the EW measure- 
ments directly; however, there is a much larger scatter on 
the high-z data. Some of this scatter will be due to a lower 
signal-to-noise ratio in the high-z data, but some might be 
intrinsic to the supernova populations that we do not ob- 
serve in the low-z sample due to its limited size. The lack 
of a well-calibrated set of supernova spectra which cover the 
blue Ca II H&K feature limits this part of the study particu- 
larly. There might be additional subtypes in the supernovae 
population, or transitional objects, which have not yet been 
identified. It is not yet clear whether there is a continual class 
of objects from underl uminous SN 1991bg-like supernovae 



(|Filippenko et al.lll992l) to ov erluminous SN1991T-like ob- 
jects (Filippcn ko et al.l Il992ll. or whether the progenitors 



are d ifferent (for example iBenetti et al . 2005; B ranch et al.l 
2006). Indeed, unusual objects objects falling outside these 
classes are now observed. 

One of the most interesting results to come out of the 
low-z trend of the Mgn feature is the fact that this feature 
seems to be constant with time and with very little scat- 
ter between objects. The feature is not made up of solely 
Mgn, but also contains Fen and Fein which are thought 
to be sensitive to temperature variations in the supernova 
and drive the relation between 7^642/443 and colour described 



in lBailev et al.l (|2009l l. Furthermore the recombination from 
Fe in to Fe n changes the opacity of the ejecta and cau ses the 
secon dary maximum in the supernova lightcurves (|Kasenl 
2006). We see weak correlations between EW Mgn and 
colour in Figure [9l but there is a large variation in the rela- 
tions between the low- and high-z data. As there is a larger 
range of colours observed in the low-z sample compared to 
the high-z sample, we would expect the any diversity to be 
seen in the low-z trend. A larger data-set at low-z would 
enable this to be investigated further. 



7 CONCLUSIONS 

The aims of this study were to observe spectroscopically 
potential Type la supernovae and confirm their type and 
redshift in order for them to be used for cosmology; to study 
the spectra of the objects themselves to verify their use as 
cosmological standard candles; and to investigate new ways 
the spectral information can be used to calibrate the objects 
for cosmology. 

The first aim is easily satisfied. Here we present obser- 
vations of 95 targets, 68 of which were confirmed as SNe la 
or probable SNe la events. 8 events were confirmed as being 
core-collapse in nature, leaving 19 observations as proba- 
bly not la in nature or of unknown type. The spectra are 
presented in Appendix [B] The spectroscopic observations of 
such a large dataset, especially when combined with other 
SNLS observations, provide an opportunity to study super- 
novae covering billions of years of cosmic look-back time. 
SNLS provides not only spectroscopic coverage, but com- 
plete light curves in four bands to fit the photometric prop- 
erties such as apparent magnitude, stretch and colour; and 
deep images of the host galaxy in five bands. In total, spec- 
tral analysis was possible for 144 objects drawn from SNLS 
survey work at the Gemini and VLT telescopes, covering a 
redshift range of 0.184 sC z s: 1.002. 

The second aim was investigated in Section 15.11 where 
high-z equivalent width measurements of the Can H&K, 
Si 11 and Mg 11 spectral features are compared to trends seen 
in a low-z sample. We find that for the Can H&K feature, 
the scatter on the measurements at low- and high-z is large, 
but the features appear to be consistent with each other. 
The Mgn feature has a much lower scatter at low redshift 
than with high-z measurements and the statistical compar- 
ison between the two sample yields Xr > L however, there 
does not appear to be a systematic shift in the mean value 
of EW Mg 11 between the two samples. The Si 11 feature has 
a lower trend at high redshift compared to low redshift. This 
is explained by the fact that E W Si 11 is seen to be anticorre- 
lated with photometric stretch. At high redshift we see more 
of the prompt component from the "A+B" model which 
are brighter and have a broader light curve due to greater 
56 Ni production. These supernovae have larger stretches and 
hence smaller EW Si 11. We also observe different correlations 
between EW Mg 11 and photometric colour for the low-z and 
high-z samples. This difference is probably due to a differ- 
ence in sample selection and the colour-cut introduced in 
the high-z target selection in particular. The origin of any 
possible relationship between colour and EW Mgn is not 
understood, but may be connected to the balance between 
Fen/Fein which are also present in the feature. It does not 
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appear to systematically affect the comparison of EW trends 
despite the high-z sample being bluer than the low-z sample. 
We find nothing within the dataset to introduce doubt as to 
the use of SNe la as standardisable candles for cosmology. 

The final aim, to exploit spectral properties to aid cos- 
mological calibration is studied in Section \5. 31 We find that 
the spectral features do not correlate with Hubble diagram 
residuals as strongly as photometric stretch or colour, but 
significant relationships are seen between the residuals and 
EW Sin and EW Mgn. This is to be expected as we see 
relations between EW Sill and EW Mgn, and stretch and 
colour, respectively. The high-z data are too noisy to be 
used to produce a spectroscopically-calibrated Hubble dia- 
gram. However, we demonstrate that such a calibration may 
be possible by defining a spectral stretch from a sub-set 
of the low-z sample based on the observed relationship be- 
tween EW Si n and photometric stretch. This is then applied 
to other low-z supernovae with known distance moduli or 
which are within the Hubble flow and we demonstrate that 
this reduces the scatter on the measurement of the absolute 
magnitude. The spectral stretch correction reduces the RMS 
on the measurement of Mb to comparable, in fact slightly 
smaller value than that obtained by photometric stretch cor- 
rection, although this is demonstrated for a small sample of 
objects. 
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Table Al. Table of observations. 1 Measured using GEMSEEING, 2 Measured using IMEXAM, 3 R150 observation, 
4 Second night of data only used 5 First night of data only used 



SN 


RA 


Dec 


UT Date 


Exposure Time 


Mode 


A cen f 


Seeing 


(SNLS name) 


(J2000) 


(J2000) 




0) 




(nm) 


(") 


05D4bo 


22:16:44.630 


-17:20:12.37 


2005-07-14 


2700 


N&S 


720 


0.43 


06Dlbt 


02:24:59:57 


-04:01:05.43 


2006-08-31 


5400 


N&S 


720 


0.45 


06Dlbz 


02:25:07.62 


-04:42:06.58 


2006-08-22 


5400 


N&S 


720 


0.52 


06Dlck 


02:26:25.37 


-04:40:43.80 


2006-08-31 


3600 


N&S 


720 


0.40 








2006-09-01 


1800 


N&S 


720 


1.03 


06Dldt 


02:27:15.20 


-04:35:44.88 


2006-09-23 


3600 


C 


680 


0.51 


06Dlgg 


02:26:28.56 


-04:44:54.29 


2006-09-23 


2400 


c 


680 


0.66 


06Dlgu 


02:24:33.91 


-04:17:37.80 


2006-11-18 


5400 


N&S 


720 


0.37 


06Dlhl 


02:26:42.95 


-04:18:22.34 


2006-12-11 


5400 


N&S 


720 


0.36 


06Dljp 


02:25:48.07 


-04:17:56.02 


2006-12-22 


7200 


N&S 


720 


0.61 


06D2iz 


10:00:34.76 


02:19:20.17 


2006-12-21 


5400 


N&S 


720 


0.64 


06D2ja 


10:00:33.51 


02:28:01.01 


2006-12-21 


3600 


N&S 


720 


0.46 


06D2ji 


09:59:51.03 


01:51:50.14 


2006-12-19 


7200 


N&S 


720 


0.66 








2006-12-20 


1800 


N&S 


720 


0.38 


06D2js 


10:02:15.70 


02:35:06.62 


2006-12-25 


7200 


N&S 


720 


1.10 2 


06D2ju 


10:01:34.65 


02:26:42.72 


2006-12-24 


5400 


N&S 


720 


0.83 


06D3ed 


14:18:47.82 


52:45:23:28 


2006-05-27 


2400 


N&S 


720 


0.82 


06D3el 


14:17:01.06 


52:13:56.84 


2006-05-27 


5400 


C 


680 


0.72 


06D3et 


14:22:04.76 


52:20:18.30 


2006-05-31 


5400 


N&S 


720 


0.56 


06D3fp 


14:20:18.62 


52:53:49.62 


2006-06-18 


2600 


C 


680 


0.55 


06D3fr 


14:19:20.24 


52:31:50.11 


2006-06-29 


3600 


c 


680 


1.07 


06D3fv 


14:21:03.31 


53:00:52.77 


2006-06-25 


1800 


N&S 


720 


0.73 








2006-06-29 


5400 


N&S 


720 


0.65 


06D3gg 


14:20:17.36 


53:10:42.27 


2007-07-20 


3600 


N&S 


720 


0.51 


06D3gh 


14:18:31.38 


52:26:34.69 


2006-06-30 


5400 


N&S 


720 


0.41 


06D3gn 


14:17:44.61 


52:21:40:42 


2006-07-06 


2700 


C 


680 


0.40 


06D3gx 


14:17:03.36 


52:56:21.82 


2007-07-27 


3600 


N&S 


720 


0.52 


06D3hd 


14:20:42.76 


52:14:20.27 


2006-07-18 


3600 


C 


680 


0.72 


06D4dh 


22:14:19.70 


-17:35:06:23 


2006-07-30 


5400 


c 


680 


0.51 


06D4dr 


22:16:46.00 


-17:20:45.49 


2006-08-19 


4800 


N&S 


720 


0.44 


06D4fc 


22:13:52.30 


-17:19:30.31 


2006-09-23 


5400 


N&S 


720 


0.39 


06D4fy 


22:14:06.22 


-17:34:19.34 


2006-09-16 


5400 


N&S 


720 


0.39 


06D4gn 


22:15:14.81 


-17:23:01.36 


2006-09-22 


7200 


N&S 


720 


0.42 


06D4hc 


22:15:40.52 


-18:09:08.71 


2006-10-13 


5400 


N&S 


720 


0.96 


07Dlcx 


02:24:58.48 


-04:01:54.23 


2007-11-10 


3600 


N&S 


720 


0.44 










3600 




715 




07Dlds 


02:25:24.44 


-04:18:49.29 


2007-11-10 


3600 


N&S 


720 


0.41 










3600 




720 




07Dldv 


02:26:51.27 


-04:47:18.29 


2007-11-10 


1800 


N&S 


720 


1.06 








2007-11-11 


1800 






0.51 










3600 




715 




07Dlca 


02:25:14:77 


-04:55:39.80 


2007-11-13 


3600 


N&S 


720 


0.45 










3600 




715 




07D2bu 


09:59:27.590 


02:33:03.48 


2007-03-15 


5400 


N&S 


720 


0.93 








2007-03-16 


1800 






0.63 








2007-03-17 


3600 






0.77 


07D2cy 


10:01:49.486 


02:40:59.09 


2007-03-22 


3600 


N&S 


720 


0.82 








2007-03-22 


5400 




715 


0.51 


07D2kc 


10:01:07.19 


02:02:28.65 


2007-12-20 


900 


C 


680 


0.67 










1800 




675 




07D2ke 


10:01:04.73 


01:59:31.81 


2007-12-20 


900 


c 


680 


0.80 










900 




675 




07D2kh 


10:00:28.94 


02:13:39.65 


2008-01-04 


3600 


N&S 


720 


0.63 










1800 




715 










2008-01-06 


1800 






0.32 


07D2ki 


10:00:26.69 


01:56:48.87 


2008-01-05 


3600 


N&S 


720 


0.54 










1800 




715 










2008-01-08 


3600 






0.54 


07D2kl 


09:59:16.95 


02:31:28.13 


2008-01-4 


3600 


N&S 


720 


0.48 










3600 




715 
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07D3ac 


14:17:33.70 


52:31:27.17 


2007-01-25 


2400 


C 


680 


0.57 


07D3af 


14:19:05.01 


53:06:08.98 


2007-01-27 


2400 


c 


680 


1.12 2 


07D3ai 


14:19:06.711 


52:40:11.69 


2007-02-14 


3600 


c 


680 


0.49 


07D3am 


14:22:36.55 


52:20:26.26 


2007-01-23 


3600 


c 


680 


0.70 


07D3ap 


14:18:10.26 


52:23:07.38 


2007-01-23 


3600 


N&S 


720 


0.54 


07D3bb 


14:18:10.31 


52:18:26.38 


2007-01-25 


7200 


N&S 


720 


0.46 


07D3bo 


14:20:58.408 


52:20:37.22 


2007-02-14 


5400 


N&S 


720 


0.54 








2007-02-15 


1800 








07D3bp 


14:21:01.32 


52:13:37.19 


2007-01-27 


7200 


B&S 


720 


0.75 


07D3bt 


14:21:25.719 


52:58:13.67 


2007-02-22 


9000 


N&S 


720 


0.55 


07D3cn 


14:16:17.255 


52:32:22.60 


2007-03-24 


3600 


N&S 


720 


0.50 








2007-03-25 


5400 




715 


0.41 


07D3cr 


14:18:12.821 


52:19:04.57 


2007-03-20 


7200 


N&S 


720 


0.41 


07D3cu 


14:20:56.531 


52:15:13.36 


2007-03-24 


2400 


C 


680 


0.49 










2400 




675 




07D3da 


14:20:06.236 


52:16:16.64 


2007-03-21 


3600 


N&S 


720 


0.36 










3600 




715 




07D3dj 


14:22:40.795 


53:02:14.69 


2007-04-12 


3600 


C 


680 


0.44 


07D3do 


14:20:35.683 


52:44:48.58 


2007-04-10 


3600 


N&S 


720 


0.54 










3600 




715 




07D3ea 


14:21:16.116 


52:40:07.88 


2007-04-12 


3600 


C 


680 


0.47 








2007-04-14 


1200 






0.66 


07D3ey 3 


14: 19:18.944 


53:04:35.92 


2007-05-14 


3600 


N&S 


720 


0.47 










3600 




710 




07D36 3 


14:21:24.800 


53:10:40.10 


2007-04-23 


3600 


N&S 


720 


0.77 










3600 




710 




07D3gm 3 


14:18:32.321 


52:43:29.62 


2007-05-19 


3600 


N&S 


720 


0.47 










3600 




710 




07D3hl 


14:16:25.147 


53:08:33.66 


2007-06-09 


3600 


N&S 


720 


0.72 










3600 




715 




07D3hu 


14:21:14.560 


52:47:53.41 


2007-06-11 


3600 


N&S 


720 


0.37 










1800 




715 




07D3hv 


14:21:05.942 


52:53:36.62 


2007-06-11 


3600 


C 


680 


0.43 


07D3hw 


14:21:48.410 


52:30:33.64 


2007-06-08 


3600 


N&S 


720 


0.53 










3600 




715 




07D3hz 


14:20:29.569 


52:23:56.79 


2007-06-10 


3600 


C 


680 


0.54 


07D3ib 


14:17:10.791 


53:04:46.91 


2007-06-16 


3600 


N&S 


720 


0.74 










1800 




715 










2007-06-18 


1800 






0.89 


07D3it 


14:18:08.730 


52:25:04.06 


2007-07-16 


3600 


N&S 


720 


0.50 










3600 




715 




07D4A 


22:16:13.16 


-18:11:46.00 


2007-11-13 


1800 


N&S 


720 


0.47 










3600 




715 




08Dlaa 4 


02:27:37.13 


-04:03:02.64 


2008-01-27 


1800 


N&S 


720 


0.69 










1800 




715 










2008-02-13 


1800 




720 


0.98 


08D2aa 5 


10:01:52.77 


02:15:57.74 


2008-01-15 


1800 


N&S 


720 


0.76 








2008-01-17 


1800 






0.95 










1800 




715 




08D2ad 


10:00:59.19 


01:45:27.31 


2008-01-09 


1800 


N&S 


720 


0.60 








2008-01-10 


3600 








08D2au 


10:02:01.84 


02:15:03.76 


2008-01-17 


1800 


N&S 


720 


0.60 










1800 




715 




08D2bj 4 


09:58:51.48 


02:07:07.38 


2008-01-10 


3600 


N&S 


720 


0.42 










3600 




715 


0.47 








2008-01-12 


1800 




715 




08D2ch 


10:01:39.61 


02:06:57.56 


2008-02-12 


1800 


N&S 


720 


0.68 










1800 




715 




08D2cl 


09:59:28.55 


02:29:49.02 


2008-02-12 


3600 


N&S 


720 


0.94 










1800 




715 




08D2dr 


10:01:54.13 


02:04:01.96 


2008-02-28 


1200 


C 


680 


0.66 










1200 




675 




08D2dz 


09:58:53.07 


01:44:19.21 


2008-03-05 


3600 


N&S 


720 


0.54 










3600 




715 




08D2eo 


10:02:03.85 


02:33:44.74 


2008-03-07 


1200 


C 


680 


0.51 










1200 




675 
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08D2fj 


10: 


00 


:57 


:55 


02: 


02 


:30. 


.68 


2008- 


■03-09 


1200 
1200 


C 


680 
675 





.61 


08D2gw 


10: 


01 


:05 


.95 


02: 


31 


:19 


93 


2008- 


-03-08 


3600 
1800 


N&S 


720 
715 


0. 


63 


08D2hw 


10: 


00 


:06 


.70 


02: 


33 


:08. 


66 


2008- 


■03-11 


1800 
3600 


N&S 


720 
715 





.69 


08D2id 


10: 


02 


:01 


.98 


02: 


13 


:13. 


38 


2008- 


-04-02 


1800 
3600 


N&S 


720 
715 





.60 


08D2iq 


09: 


58 


:56 


.05 


02: 


:34 


:26. 


.99 


2008- 


-04-08 


1800 
3600 


N&S 


720 
715 





.45 


08D2kj 


09: 


59 


:13 


.37 


02: 


■25 


:04. 


45 


2008- 


-04-08 


3600 
3600 


N&S 


720 
715 





.43 


08D3au 


14: 


19 


:43. 


.14 


53: 


06 


:31, 


.22 


2008- 
2008- 


-02-15 
■02-17 


1800 
3600 


N&S 


720 
715 






.61 
.51 


08D3bh 


14: 


17 


:06 


.50 


53: 


01 


:37 


.38 


2008- 


-02-15 


1800 
1800 


N&S 


720 
715 





.69 


08D3dc 


14: 


20 


:44. 


.79 


52: 


46 


:48, 


.15 


2008- 
2008- 


-03-08 
-03-10 


3600 
3600 


N&S 


720 
715 






.51 
.55 


08D3dx 


14: 


17 


:47. 


.82 


52: 


■23 


:02 


76 


2008- 


-03-08 


3600 
3600 


N&S 


720 
715 





.58 


08D3fp 


14: 


16 


:28. 


.53 


53: 


06 


:10 


23 


2008- 


■04-08 


3600 
3600 


N&S 


720 
715 





.64 


08D3fu 


14: 


22 


:35 


.86 


53: 


00 


:36 


98 


2008- 


-04-02 


1800 
1800 


N&S 


720 
715 





52 


08D3gb 


14: 


21 


:53 


.41 


52: 


11 


:47, 


89 


2008- 


-04-03 


1200 
1200 


C 


680 
675 





.67 


08D3gf 
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Table A2. Classifications and redshifts. Objects types arc defined in Section 12.31 
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Target Type Redshift Redshift Error CI Redshift Source 
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3 
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5 
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5 
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o 


host 
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SNIa 


0.451 
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5 


host 
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2 




07D3bo 


SNIa? 


0.92 


0.02 


3 


target 


07D3bp 


SN? 


0.769 


0.001 


2 


host 


07D3bt 


SNIa? 


0.91 


0.01 


3 


target 


07D3cn 


SNIa 


0.898 


0.001 


4 


host 


07D3cr 


SNIa 


0.746 


0.001 


4 


host 


07D3cu 


SNIa 


0.512 


0.001 


5 


host 


07D3da 


SNIa 


0.837 


0.001 


1 


host 


07D3dj 


SNIa 
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0.001 


5 


host 


07D3do 


SNIa? 
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0.01 


3 


target 


07D3ea 


SNIa 


0.471 


0.001 


5 


host 


07D3ey 


SNIa 
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4 


host 


07D3fi 
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07D3gm 


SNIa 


0.83 


0.02 


5 
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SNIa 


0.67 
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5 
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Table A2 - continued 



07D3hu 


SNIa 


0.572 


0.001 


4 


host 


07D3hv 


SNIa 


0.351 


0.001 


5 


host 


07D3hw 


SNIa 


0.748 


0.001 


5 


host 


07D3hz 


SNIa 


0.506 


0.001 


5 


host 


07D3ib 


SNIa 


0.681 


0.001 


4 


host 


07D3it 


SNIa 
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0.001 


4 


host 


07D4A 


SNIa 


0.503 


0.001 


5 


host 


08Dlaa 


SNIa? 


0.54 


0.01 


3 


target 


08D2aa 


SNIa 
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5 


host 


08D2ad 


SNIa 


0.554 
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5 


host 


08D2au 
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2 


host 
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SNIa 
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4 
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SNIa? 
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3 


host 
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SNIa? 
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0.001 


3 


host 


08D2dr 


SNIa 


0.355 


0.001 


4 


host 


08D2dz 


SNIa 


0.65 


0.01 


4 


target 
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SNII 


0.125 


0.001 





host 
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0.001 


1 


host 
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SNIa 
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0.001 


5 


host 
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SNIa 


0.746 


0.001 


4 


host 
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0.001 


4 


host 
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0.001 


5 


host 
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host 
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5 
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5 
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5 
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SNIa 
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5 
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Table A3. EW results for objects presented in this paper. 
Objects marked * arc removed from further analysis. 



Object 
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(A) 


(A) 


DR'D'Wl 






a 1 Q-Ufi 1 R 


yu.o 1 nz io.uz 
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1U1. r It 1 1 . 


V7l.)J_-'t>i_- L 


89 
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89 
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61±6.49 
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62 
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79 
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31.9±2.06 
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07D3hv 
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93±4.96 
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70.7±1.05 
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07D3hz 
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10.96±1.60 
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07D3ib 
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62±43.12 


17.9±21.22 
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07D3it 


49. 


31±14.42 


28.04±9.76 


147.35±16.35 
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Table A4. EW results for the Bronder sample. Objects 
marked * are removed from further analysis 
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1 ns S4+1 fi 28 


21 1 Q+1 3 RQ 

Z -L . A- y _l LJ.Uc; 
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o'-t. -Luznu . 00 
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lUiitUl 1_ J_ rt • iiU 


n5m hv 
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oy . -Lo_i_ y .iiu 


05Ti1 rr 


40 Afi-Ui 1 on 
10 .uuzrz ii.yu 


1 7 57-U5 77 


QI 4C-L1 4 Q7 
.L .iozc ii.y / 


05m no 


54 Q7-U1 1 Q5 


n f|Q+4 fiQ 
u. uyzci. uy 


1 04 43+1 9 99 

-LUrL-.rtOZIZ _LZi. ii 


05! VI OTTI 


1 99 48+30 87 


99 3Q_L99 07 


1 97 09+57 98 

-LZi / .UZrZLZO 1 .^O 


05r>9ah 




11 qi -Li n 98 

_L J- .O -L TJZ -LU . ZiO 


7Q 39-UI 9 1 fi 
1 y .ozzn j-Zr. iu 


\J > J 1_J Z (XI 1 






QI fiO+1 8 fiO 

^ J. . ■ H J .LO. UU 




«q 55-UI 4 50 


94 40+fi 89 


1 48 44+1 3 1 fi 

14:0- H-'H- 1 J_0. -LU 


05D2ia 

U ' J 1— » iij CI. 




« qfi+'i S7 

U.OU_LO.O 1 


1 q« c:q+a 1 2 

-LOU. Oy_LO . -Lii 


UU-L-'ZllL 


1 n5 5+14 77 


4 3Q-U5 04 
1. oozco. yi 


1 94 1 9_|_i 3 QQ 
±z i ±. j-Zi— L7. .Lo. oy 


W ' J 1 ' i. V 7 


1 01 98+1 Q 1 


(S 09+1 9 83 

U. L/ii_l Li.OJ 


121 q2+31 37 

lii. y ii— i— i.ji 


UO-L-'OCLA. 


1 1 Q Q7_|_Q 7fi 


5 f|7+4 84 

O. U 1 ZUrt. Ol 


75 1 7+8 3fi 
1 j.i 1 zno .ou 


05r»3rf 




1 q 79-Ui 9 1 7 
-Lo . / zzn iz. 1 1 


1 49 31 _|_i 3 «3 

lii. O -LZ L7_ .LO. UO 


UU-L-'OLl 


70, Q7_L99 47 


1 1 0Q+5 fi5 
_L j. . uyzco.uo 


71 99-UQ fi4 


05D3cx 


56.10±13.18 


18.24±10.13 


110.71+20.46 


05D3jq 


55.17±16.18 


8.61±5.30 


18.52+5.63 


05D3kt 


127.22±13.46 


13.05±7.17 


83.19+11.43 


05D31b 


115.57±10.36 


9.25±4.47 


88.22+9.21 


05D3mh 


74.90±31.00 


17.50±9.68 


147.49+23.48 


05D3mq 




36.46±34.80 


117.68+16.37 


05D3mx 


45.99±22.97 


19.45±8.79 


108.38+11.20 


05D4bm 


144.93±7.48 


12.60±3.43 


75.91+4.910 


05D4dt* 




2.13±2.64 


6.38+3.17 


05D4dy 


95.96±25.01 


23.12±11.81 


125.51+23.31 


05D4fo 




28.59±4.52 


80.82+5.90 


05D4gw 


143.8±27.02 


12.85±10.66 


131.23+21.20 


06D3bz 


113.9±21.40 


13.48±9.22 


81.48+24.17 


06D3cn 




6.57±2.6 


103.29+0.93 
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Table A5. EW results for the VLT sample. Objects marked 
* arc removed from further analysis. 



Object 


Can 


Si II 


Mgn 




(A) 


(A) 


(A) 


UO.L. / _1_ Ul 
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O 1 .OOIOU.H:U 


i 4 Q0-I-1 t 49 


1 11 44-T-34 S 1 ^ 

IOO . 11IOH:.OJ 


n^m m 


i nQ nQ-r-^s n7 


9 97c_li ft «4 
Zr. z i oxc iu.U4: 


1 9c: QK_|_4n no 


n^m Ht 

UO-L' 1111 


«7 QCT_L1 4 Q1 


lO.OUX 1 .OO 


1 11. zoic 10.00 


03m fc 


1 fi4 Qo+o 4Q 

IUt:. yU_ l_ O . iy 


O Q4-L-Q nn 

O . ■ n 1 ■ > . \ J\J 


i ni ',7-l-fi 98 


03D1A 


QA Q1 _|_1 1 1 Q 

y i . o i nz i i . i o 


U.U4:X4:.UO 


04 QQ+i 9 99 

OL ■ I' * 1 iz. zz 


UO-L-'idg 


7Q .1 +9 41 


4.61±1.43 


«q 74+q nfi 

uy. 1 i_i_o.uu 


03f)4at 

UO-L-'idL 


•J 1 . U J — 1 LO . UU 


cc (S7+Q 1 9 

J.UI _1_ J. li 


qn 27+28 34 




9,7 ni _L1 K Q1 
Ol .Ul ZL -LO .0 1 


1 4 97-uq nn 

i4:.zi xcy.ou 


Q1 70-1-1 7 4^ 
01.1 um 1 1 •'-to 


UOl^ILl y 


1 nfi ^n-n r ia 


o.O'Jiy.u i 


07 nQ_|_9n 1 4 
y i .uorxizu. 11 


UH-* J_U.(_- 


J. .LO . yorxrz.iy 


9fi 7f|-r-1 t^ 
ZU. 1 OXC 1 . OO 


fit; 9«-|_9 40 


04DlfT 


1 R7 7 c i-r-^9 ^fi 
iui . / orxroz.ou 


4-O.U'iX/nJ.UU 




D4D1 hx 


liU . OO— 1 LU.UJ 


91 nc;_L i n on 

— 1 .OO—I LU.OU 


1 04 1 9+1 1 So 
iui. 1 ^ 1 1 1 . 00 


UtJJ 11 V 


79 79+44 1 n 

1 Z. 1 ZrXT1:1:.lU 






04Dlkj 


1 07 7fi+X Q1 

IU 1 • 1 U_ l_ O. y 1 


8 22+4 99 


88 q8+1 20 

00. yo_i_ iu.iW 




1 r:n 7Q-L1 9 on 

Lijy. / O^ LZ.OU 




1 QQ CCQ_|_9Q «Q 

100. oyrxizo. uo 


D4D1 nw 

UlL/ 1U W 


1 OQ Q4+1 1 71 

J_ C/ . £/rr_l Lit t -L 


.9 ni+1 1 fS4 






7fi Q4_LQ9 79 


ZrO . yoxcoz .ou 




n4m np 


UK) . OO—I LL/.OO 


« 98+4 71 


on qc;+i n 21 
oy . y u_i_ iu.ii 


UH 1 1 ]->]J 


Q7 OQ-U1 9 91 


1 k 1 Q-Ufi c;9 

10 . IOXCU. oz 


09 1 9_|_i x. crn 
oz. 1 zzc 1 0. ou 


04ni nd 


on 23+1 fi 87 


oq ro_li n 4c; 




n4ni rh 

UH-* 1L 11 
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Table A6. Photometric properties of the low-z sample fitted using SiFTO 



SN 


Redshift 1 


Stretch 


Colour 


B-band 2 
(mag) 


Independent Distance 


SN 1 QSQR 

DIN 


0.004 


n qn^-i-n 097 


0.374±0 


.029 




QO ^0 -1- 0Q a 


din j_yyuiN 


o ooa 


i 009-1-0 01 


0.051±0 


.013 


1 o f!Q-l_o 09 


39 90 -1- OQ a 


SN 1 QQ4T) 


0.003 


n ai 9-1-0 OOzL 


-0.104±0 


.006 


1 1 7^-1-0 09 

J- J- . / OZL\KJ . UZ 


91 08 -1- 90 b 


SN 1 QQfiX 
din lyyuA 


0.008 


U.OU'lxU.Ul / 


-0.023±0 


.016 


1 9 QQ-l-0 04 




OJ.N _Lfi70cll-£ 


0.004 


n Qc:e;_|_n nnq 


-0.146±0 


.009 


1 9 *31 -)_o 01 


Ol.O/ ZIZ U.J-O 


SN 1 QQShn 

DIN Ic/cJOULl 


0.004 


n qc;4_|_n ni « 
u. yj'ixu.u lu 


0.262±0 


Oil 


1911 _)-0 09 

A-Zt . A. L ZIZ \J . Viz 


90 94 -1- 1 1 a 


SN 1 QQQpp 


0.011 


1 OR^-r-O OOfi 


0.261±0 


.006 


1 A 8^-1-0 01 




SN 9000F 


0.004 


i OfiD-l-n 01 9 

J. .UUUXU.U1 ^ 


0.128±0 


.035 


1 9 ac;-|-o 1 ^ 




SN 9009hn 


0.005 


n Q40-|-n 01 1 


0.379±0 


.014 


1 q qc^-l-O 09 
lo.yjiu. viz 




SN 9009nr 
DIN ZUUZt!! 


00Q 


Q7fi-Un 00Q 


0.124±0 


.019 


1 A 9^-1-0 Ofi 




SN 2003du 


0.007 


1.006±0.007 


-0.112±0 


.008 


13.46±0.01 


_ 


SN 2004eo 


0.015 


0.869±0.005 


0.048±0. 


.013 


15.07±0.05 




SN 2004S 


0.010 


0.976±0.034 


0.011±0 


.025 


14.14±0.05 




SN 1997do 


0.010 


0.944±0.029 


0.002±0 


.027 


14.27±0.05 




SN 1998V 


0.017 


1.001±0.039 


-0.027±0. 


.025 


15.07±0.08 




SN 1998dh 


0.008 


0.907±0.017 


0.029±0 


.025 


13.85±0.04 




SN 1998ec 


0.020 


0.975±0.070 


0.115±0 


.058 


16.16±0.12 




SN 1998eg 


0.024 


0.919±0.065 


-0.010±0 


.022 


16.11±0.05 




SN 1999cc 


0.032 


0.784±0.016 


0.005±0 


.010 


16.78±0.02 




SN 1999gd 


0.019 


0.920±0.056 


0.365±0. 


.024 


16.94±0.04 




SN 1999gp 


0.026 


1.167±0.015 


0.031±0 


.012 


16.04±0.03 




SN 2000fa 


0.022 


1.027±0.031 


0.043±0 


.027 


15.85±0.05 




SN 2001V 


0.016 


1.128±0.012 


0.050±0 


.015 


14.62±0.02 





Relative to CMB, 2 At date of maximum. 



References: a ISaha et al] <2006fl 6 lAihar et all ll200ll) 
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Table A7. EW measurements for the low-z spectra. 



Object 


Phase 


Can 


Sin 


Mgn 




(Days) 


(A) 


(A) 


(A) 


SN 1989B 







26.56 


98.58 


SN 1989B 


6 


— 


39.82 


128.16 


SN 1989B 


11 


— 


— 


172.61 


SN 1990N 


-8 


— 


6.92 


92.35 


SN 1990N 


2 


- 


15.44 


80.78 


SN 1994D 


-11 


193.27 


10.09 


59.06 


SN 1994D 


-9 


— 


15.13 


68.78 


SN 1994D 


-8 


— 


20.16 


57.26 


SN 1994D 


-5 


100.77 


20.10 


65.99 


SN 1994D 


-4 


96.37 


18.09 


74.21 


SN 1994D 


-2 


99.63 


23.05 


79.46 


SN 1994D 


2 


63.20 


23.45 


71.85 


SN 1994D 


4 


75.34 


21.47 


65.94 


SN 1994D 


10 


— 


20.39 


90.39 


SN 1994D 


11 


65.92 


16.71 


108.98 


SN 1996X 


-2 


41.29 


17.18 


86.65 


SN 1996X 





81.01 


17.77 


92.88 


SN 1996X 


1 


82.05 


17.75 


91.66 


SN 1996X 


7 


80.56 


13.71 


93.10 


SN 1998aq 


-9 


29.61 


10.29 


92.22 


SN 1998aq 


-8 


58.20 


10.27 


95.58 


SN 1998aq 


1 


38.40 


10.18 


79.06 


SN 1998aq 


2 


38.79 


12.37 


81.56 


SN 1998aq 


3 


37.22 


12.01 


94.22 


SN 1998aq 


4 


38.15 


12.78 


96.98 


SN 1998aq 


5 


40.50 


12.70 


94.88 


SN 1998aq 


6 


40.15 


12.61 


96.50 


SN 1998aq 


7 


39.76 


12.88 


95.41 


SN 1998bu 


-4 


91.48 


15.24 


99.23 


SN 1998bu 


-3 


89.16 


15.24 


99.23 


SN 1998bu 


-2 


87.53 


14.30 


98.41 


SN 1998bu 


8 


81.97 


15.27 


90.43 


SN 1998bu 


9 


83.49 


16.36 


91.26 


SN 1998bu 


10 


86.99 


14.78 


96.61 


SN 1998bu 


11 


82.38 


13.68 


112.96 


SN 1999ee 


-11 


95.76 


1.14 


106.85 


SN 1999ee 


-9 


98.16 


1.75 


93.78 


SN 1999ee 


-4 


103.10 


7.29 


104.13 


SN 1999ee 


-2 


108.93 


7.92 


102.26 


SN 1999ee 





109.82 


8.91 


105.23 


SN 1999ee 


5 


129.76 


18.67 


112.07 


SN 1999ee 


7 


— 


20.11 


114.87 


SN 1999ee 


9 


132.32 


13.38 


114.11 


SN 2000E 


-9 


148.74 


9.88 


105.79 


SN 2000E 


-5 


— 


9.97 


82.18 


SN 2000E 


5 


142.57 


12.86 


91.46 


SN 2002bo 


-11 


195.34 


20.38 


107.06 


SN 2002bo 


-6 


152.29 


17.78 


105.13 


SN 2002bo 


-5 


191.09 


23.85 


109.45 


SN 2002bo 


-4 


137.57 


22.08 


104.52 


SN 2002bo 


-3 


137.96 


22.85 


104.29 


SN 2002bo 


-2 


142.47 


33.86 


114.94 


SN 2002bo 


-1 


145.12 


23.34 


108.19 


SN 2002bo 


4 


118.94 


36.97 


117.07 


SN 2002er 


-11 


170.87 


23.61 


107.01 


SN 2002er 


-9 




27.09 


102.27 


SN 2002er 


-8 


168.78 


22.12 


98.25 


SN 2002er 


-7 


179.28 


39.70 


99.45 


SN 2002er 


-6 


161.45 


20.32 


102.01 


SN 2002er 


-5 


162.72 


20.72 


103.20 


SN 2002er 


-4 


171.20 


23.61 


105.90 


SN 2002er 


-3 


152.98 


19.58 


91.61 


SN 2002er 


-2 


152.98 


19.58 


91.61 
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Table A7 - continued 



SN 2002cr 


-1 


149.57 


22.05 


97.47 


SN 2002er 





- 


21.51 


115.32 


SN 2002cr 


2 


- 


23.78 


97.44 


SN 2002er 


4 


111.98 


25.90 


99.99 


SN 2002er 


5 


- 


35.25 


128.83 


SN 2002er 


6 


108.88 


24.68 


88.57 


SN 2002cr 


10 


108.03 


8.29 


110.61 


SN 2003du 


-11 


145.34 


17.62 


82.77 


SN 2003du 


-8 


176.57 


14.33 


92.39 


SN 2003du 


-7 


155.24 


12.39 


93.97 


SN 2003du 


-6 


134.42 


10.66 


82.10 


SN 2003du 


-4 


133.97 


12.75 


83.85 


SN 2003du 


-2 


132.71 


12.15 


89.28 


SN 2003du 


-1 


138.70 


12.40 


90.29 


SN 2003du 





136.87 


13.52 


93.04 


SN 2003du 


1 


141.85 


14.25 


93.59 


SN 2003du 


3 


134.86 


14.39 


80.77 


SN 2003du 


4 


119.98 


14.63 


79.39 


SN 2003du 


7 


89.99 


16.42 


85.42 


SN 2003du 


8 


89.31 


15.90 


82.90 


SN 2003du 


9 


86.63 


17.18 


84.06 


SN 2003du 


10 


85.06 


16.99 


88.40 


SN 2004dt 


-10 


81.19 


43.41 


118.9 


SN 2004dt 


-9 


89.76 


40.28 


119.82 


SN 2004dt 


-7 


91.64 


34.27 


109.54 


SN 2004dt 


-6 


97.04 


30.46 


115.12 


SN 2004dt 


-4 


97.79 


27.33 


113.16 


SN 2004dt 


-3 


96.84 


32.15 


130.73 


SN 2004dt 


-2 


98.46 


24.45 


120.35 


SN 2004dt 


-1 


116.19 


20.96 


104.04 


SN 2004dt 


2 


120.69 


25.87 


99.30 


SN 2004dt 


3 


131.45 


24.94 


115.57 


SN 2004co 


-6 


136.35 


24.30 


107.61 


SN 2004co 


-3 


118.77 


32.41 


128.00 


SN 2004eo 


2 


96.18 


25.44 


104.67 


SN 2004co 


7 


92.51 


15.36 


130.34 


SN 2004S 


1 


148.99 


10.46 


110.43 


SN 2004S 


7 


144.40 


13.59 


121.01 


SN 1997do 


-11 


- 


14.33 


105.19 


SN 1997do 


-10 


- 


12.12 


97.86 


SN 1997do 


-7 


- 


20.44 


113.20 


SN 1997do 


-6 


- 


15.18 


103.68 


SN 1997do 


8.5 


- 


11.52 


101.07 


SN 1997do 


10.5 


- 


11.94 


123.93 


SN 1998V 


0.5 


- 


14.54 


92.82 


SN 1998V 


1.5 


- 


16.38 


103.11 


SN 1998V 


2.5 


- 


13.64 


90.24 


SN 1998dh 


-9 


- 


24.87 


102.04 


SN 1998dh 


-8 


- 


26.35 


105.01 


SN 1998dh 


-7 


- 


24.21 


105.78 


SN 1998dh 


-5 


- 


27.31 


110.46 


SN 1998dh 


-3 




30.62 


107.46 


SN 1998dh 





- 


27.46 


109.28 


SN 1998cc 


-2.5 




23.23 


115.56 


SN 1998cc 


-1.5 




18.02 


109.87 


SN 1998eg 







22.23 


98.36 


SN 1998eg 


5 




24.91 


106.30 


SN 1999cc 


-3 




32.39 


107.29 


SN 1999cc 


-1 




32.82 


112.07 


SN 1999cc 


0.5 




26.97 


116.27 


SN 1999cc 


2 




34.46 


108.90 


SN 1999ej 


-0.5 




28.54 


91.06 


SN 1999ej 


2.5 




27.16 


91.83 


SN 1999ej 


4.5 




22.25 


90.77 


SN 1999ej 


9 




16.85 


123.35 
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Table A7 - continued 



SN 1999gd 


2.5 


36.35 


117.13 


SN 1999gd 


9.5 


26.53 


104.11 


SN 1999gp 


-4.5 


6.28 


111.09 


SN 1999gp 


-1.5 


5.55 


107.13 


SN 1999gp 


0.5 


10.44 


89.89 


SN 1999gp 


3 


9.64 


109.73 


SN 1999gp 


5 


22.98 


125.49 


SN 1999gp 


7 


24.69 


114.61 


SN 2000fa 


-10 


12.55 


102.34 


SN 2000fa 


1.5 


25.13 


120.30 


SJN iOOOia 


2.5 


17.96 


100.17 


SN 2000fa 


4.5 


19.19 


105.09 


SN 2000fa 


9.5 


25.16 


99.18 


SN 2001V 


-11 


10.01 


82.89 


SN 2001V 


-10 


7.81 


79.05 


SN 2001V 


-9.5 


7.98 


81.46 


SN 2001V 


-7.5 


12.15 


96.19 


SN 2001V 


-6.5 


12.07 


98.80 


SN 2001V 


-5.5 


10.11 


94.55 


SN 2001V 


-3.5 


9.10 


89.31 


SN 2001V 


9.5 


15.82 


92.23 


SN 2001V 


10.5 


19.66 


93.89 
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APPENDIX B: HIGH-Z SPECTRA 

Here we present all the SNLS supernova candidates observed 
at Gemini between May 2006 - May 2008. The spectra show 
the raw data in light blue with a smoothed version in dark 
blue. The smoothing is done using a Savitzky-Golay smooth- 
ing filter (jSavitzkv fc Golavlll964l ) . The spectra are grouped 
according to their classifications. In the data plots, if any 
separate host galaxy extraction was possible it is plotted in 
green. Any narrow lines from host emission/absorption are 
marked. The results of galaxy template subtraction from the 
method presented in Section \3. 31 are plotted in red. 
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Figure Bl. Type la Supernovae (CI = 4 or 5) 
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Figure B2. Type la Supernovae (CI = 4 or 5) 
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Figure B3. Type la Supernovae (CI = 4 or 5) 
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Figure B4. Type la Supernovae (CI = 4 or 5) 
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Figure B5. Type la Supernovae (CI = 4 or 5) 
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Figure B6. Typo la Supernovae (CI = 4 or 5) 
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Figure B7. Type la Supernovae (CI = 4 or 5) 
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Figure B8. Probable Type la Supernovae (CI = 3) 
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Figure B9. Probable Type la Supcrnovae (CI = 3) 
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Figure BIO. Unknown (CI=2) 
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Figure Bll. Unknown (CI=2) 
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Figure B12. Probably or Definitely Not Type la Supernovae (CI=1 or 0) 
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Figure B13. Probably or Definitely Not Type la Supernovae (CI=1 or 0) 



